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N , SUMMARY 

-- ..-4 

The computer program for a comprehensive analytical model 
of rotorcraft aerodynamics and dynamics is described. This 
analysis is designed to calculate rotor performance, loads, and 
noise; the helicopter vibration and gust response; the flight 
dynamics and handling qualities; and the system aeroelastic stability. 
The analysis is a combination of structural, inertial, and 
aerodynamic models, that is applicable to a wide range of problems 
and a wide class of vehicles. The analysis is intended for use in 
the design, testing, and evaluation of rotors and rotorcraft, and 
to be a basis for further development of rotary wing theories. 

This report documents the computer program that implements the 
analysis . 
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The total length of the block is given in parentheses after the label. 
Then all variables in the block axe listed. The left-hand column gives 
the variable name, and the right-hand column gives the location of the 
variable in the common block. Finally, a description of the variable is 
provided (except for variables in blocks with labels of the form xxDATA, 
which are input parameters). Only the common blocks for rotor are 
described; the common blocks for rotor ff2. have an identical structure. 
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RTR1CM(1070) 


OMEGA 

rotor speed fl, (rad/sec) 

1 

MTIP 

tip Mach number flR/cg 

2 

GAMMA 

Lock number ^ 

3 

CMEAN 

mean chord Cgj 

4 

IB 

characteristic inertia 

5 

NBM 

number of bending modes 

6 

NTM 

number of torsion modes 

7 

NGM 

zero if no gimbal or teeter mode 

8 

NBHT 

number of mean bending deflection modes 

9 

GLAG 

®lag n 

10 

MLD 

11 

DZLD 


12 

CGC 

^GC " S’ “ 

13 

CGS 


14 

NUGC 

■^C 

15 

NUGS 

■^GS 

16 

CTO 

collective control damping C^/l^n 

17 

CTC 

cyclic control damping C^ A-b’^ 

18 

CTO 

rotating control damping 

19 

RA(30) 

aerodynamic radial stations, r^^, i = 1 to MRA 

20 

ERA(30) 

aerodynamic segment length Zlkr^^, 1 = 1 to MRA 

50 

FTIP(30) 

tip loss multiplicative factor f. , 1 •= 1 to MRA 80 

PSI21M 

(rad), 0. for rotor #1 (for BODYM, 
M0TNH7 WAKEN, ENGNM) 

110 

PSI21W 

. (rad), — for rotor #2 (for 

waken; wakeg) 

111 

MUX 

A 

112 

MUY 


113 

MUZ 

Az 

114 

RGUST(3.3) 

"^G 

115 

CHUB(6,16) 

c 

124 

cbhub(3,3) 

c" (including factor 

220 

CHUBT(16,6) 

T 

c 

-24- 

229 
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ALFHP 

(deg) 

RTRIGM 

325 

PSIHP 

*^HP 

326 

MAT 

M + 

at 

327 

CD( 2 ) 

C- for drive train 
D n 

328 

CPSI( 2 ) 

for drive train motion 

330 

pinter( 36 ) 

burst tip vortex in wake model 

“^inter ^3* J = 1 to MPSI 

332 

pburst( 36 ) 

(rad) at Vj. 0 = 1 to MPSI 

368 

eixxb( 5 i) 

inertial and structural data at 
r = e + (j-l)Ar, j = 1 to MRB +1 

XX 

/ 4-04 

EIZZB( 51 ) 

rL^^/Ei 
' zz 

455 

MASSB( 51 ) 

m 

506 

’IVISTB( 51 ) 


557 

CENT( 51 ) 


608 

ITHETB( 51 ) 

inertial and structiiral data at 
r = r^^ + ( j-l)^r, j ■= 1 to MRB +1 

\ IL 

659 

gjb( 51 ) 

xt5rVgj 

710 

MASSl( 5 l) 

inertial data at 
r = (j-l)^r, j = 1 to MRM +1 

mnVi^ 

761 

ITHETl( 5 l) 


812 

XII( 51 ) 

Xj/R 

863 

XCI( 51 ) 

x^/R 

914 

TWISTl( 5 l) 

% 

kpVB^ 

965 

KP 2 I( 51 ) 

1016 

IPITCH 

2 2 

blade pitch inertia (slug-ft or kg-m ) 

1067 

KTO 

control system stiffness 
(ft-lb/rad or m-N/rad) 
collective 

1068 

KTC 

cyclic 

1069 

KTR 

reactionless 

1070 
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RHlCMf 12792) 

hrtr(i 6, 16,21 ) 
hbody(i6,6,io) 

HENG(6,10) 


complex rotor transfer function matrix, j 
size NBM+NTM+NGM; n = 0 to HHARM " 

complex airframe transfer fvuiction matrix, 

H'^c^ ; n = p = 1 to 

MHARMF 

complex drive train transfer function matrix, 

”n ^ ! n = pNnAlj.gf . P = 1 to 

MHARMF 


1 

10753 

12673 
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B0rBfCM(446) 

AM0I®1(6,16) 

C^k» rotor #1 hub (dimensionless) 

1 

AM0I®2(6,16) 

^Tk’ rotor #2 hub (dimensionless) 

61 

KMSTCl(lO) 

pitch/mast-bending coupling (dimensionless) 
Kj^Cj^ for rotor #1 

121 

KMSTSl(lO) 

^0^ rotor #1 

131 

KMSTC2(10) 

%Ck rotor #2 

141 

KMSTS2(10) 

KjjSjj for rotor #2 

151 

ADAMPA(IO) 

aerodynamic damping (2'^/v-aA)(q/v)Fqk^ 

l6l 

ACNTRL(4,10) 

control derivatives (2^ /«~ 2 Jl)qFq^^ 

171 

AMASS( 10) 

* 

211 

ADAMPS(IO) 

\*«s^k 

221 

ASPRNG(IO) 

K ^k 

231 

MSTAR 

M* 

241 

MSTARG 

M*g 

242 

istar(3,3) 

I* 

243 

CWS 

Cy/v- = (a/2 2() M*g 

252 

HMASS 

aircraft mass (slug or kg) 

253 

HAM 

number of airframe modes 

254 

RSF10(3,3) 

aircraft description ( = O) 

R„p for rotor #1 

or 

255 

RSF20(3,3) 

R_„ for rotor #2 
br 

264 

RHUB10(3) 

"r at rotor #1 hub 

273 

RHUB20(3) 

at rotor #2 hub 

276 

RWB0(3) 

T at wing/body 

279 

RHT0(3) 

"r at horizontal tail 

282 

RVT0(3) 

at vertical tail 

285 

ROFFO(3) 

"r off rotor 

288 
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BODYCM 


RSF1(3,3) 

aircraft description 
Rgp for rotor #1 

291 

RSF2(3»3) 

Rgp for rotor #2 

300 

RHUB1(3) 

r at rotor #1 hub 

309 

RHUB2(3) 

”r at rotor #2 hub 

312 

RWB(3) 

at wlng/body 

315 

RHT(3) 

r at horizontal tall 

318 

rvt(3) 

r at vertical tail 

321 

ROFF(3) 

off rotor 

324 

TCFE(11,5) 

T 

CFE 

327 

VXREXF(3) 

(V X ) R^ kp 

382 

mvxre(3,3) 

— M* (Vx ) Rg 

385 

GMTRX(3,3) 

G 

394 

IBOI]Y(3,3) 

R^I* R 
e e 

403 

REULER(3,3) 

R 

e 

412 

RFV(3.3) 

^FV 

421 

RFE(3.3) 


430 

KE(3) 


439 

VELF(3) 

V 

442 

VCLIMB 

^ climb 

445 

VSIDE 

V XJ 

side 

446 
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ENGNCH(131) 


QTHRTL 

IBNG 

KMIl 

KMI2 

KMR 

KMEl 

KME2 

KPGOVE 

KPGOVl 

KR30V2 


^E 

^MIo 


^MR 

*^ME2 


governor proportional gains , 
engine 
rotor #1 
rotor #2 


Kp*r2- 


NDH 

TIGOVE 

TlGOVl 

T1GOV2 

T2GOVE 

T2GOV1 

T2GOV2 

QEDAMP 

IRSTAR 

MENG(6,6) 

SEHG(6,6) 

DENG(6,6) 

HBNG0(2,2) 


number of drive train modes 

governor time lag# 
engine 
rotor #1 
rotor #2 

governor time lag, 
engine 
rotor #1 
rotor #2 

O * 

I* 

•^R 


mass matrix for H, 


-1 


-1 


n 


-1 


spring matrix for H 

damping matrix for H 

-1 ‘ 

Hq for static elastic motion 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
56 
92 

128 
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GUSTCMC 12989) 


gxist components, velocity axes 


VGWBV(3) 

at wing/t>ody, ^ 

1 

VGHTV(3) 

at horizontal tail, 

4 

VGVTV(3) 

at vertical tall, g^ 

7 

vGR1V(3,30,36) 

at rotor #1, gCr^^, 

10 

vgR2v(3,30,36) 

at rotor #2, g(r,, y.) 

3250 

VGHUB1(3) 

at rotor #1 hub, g (for waice geometry) 

6490 

VGHUB2(3) 

at rotor #2 hub, "g (for uake geometry) 
gust components, F axes 

6493 

VGWBF(3) 

at wing/body, ^ 

6496 

VGHTF(3) 

at horizontal tail, g^ 

6499 

vgvtf(3) 

at vertical tail, g^ 
gust components, S euces 

6502 

vgris(3,30,36) 

at rotor #1, ^r. , V-) 

6505 

VGR2S(3,30,36) 

at rotor #2, g(r, , s>.) 

^ 3 

transient control 

9745 

VPraAN(5) 


12985 
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C0NTCM(32) 

VCNTRL(ll) 

THETFT 

PHIFT 

THETFP 

PSIFP 

THETAT 

PSIT 

DVB0EY(6) 

DOMEGA 

DDZF 

VPIL0T(5) 

TGOVRl 

TGOVR2 


control vector (rad): ■p 

"v = (©7g ©I* ©-rs ie ®t) 

rotor#l rotor#2 airframe 


(rad) 

12 

‘H’T 

13 

e^p (rad) 

14 

Vpp (rad) 

15 

©.p (rad) 

16 

% (rad) 

17 

airframe motion (dimensionless) 

18 

% ‘^f ^ yF 

(static; dimensionless) 

s 

24 

Zp (dimensionless) 

25 

pilot control vector (rad); 

26 

""p = ( ^ ^ ^ ^ ^t^ 

^^^govr^rotor#! (^*^^ 

31 

^^^govr^rotor#2 (^^) 

32 
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C0NVCM(80) 


mean squzire motion (rotor #l) 


BIMS(IO) 


1 

T1HS(5) 

e 

11 

BGIMS 


16 

P1MS(16) 


17 

PS1MS(6) 

mean sqxiare motion (rotor #2) 

33 

B2MS(10) 

P 

39 

T2MS(5) 

B 


BG2MS 


54 

P2MS(16) 

4> 

55 

PS2MS(6) 

T 

71 

GIMS 

mean square circulation (rotor #l) 

77 

G2MS 

mean square circulation (rotor #2) 

78 

COUNTM 

integer peirameter: number of motion iterations 

79 

COUNTC 

Integer parameter t moaber of circulation 
iterations 

80 



MD1GM(6773) 


T750LD 

old (initialized to 1000.) 

1 

NBMOLD 

old NBM (initialized to O) 

2 

NTMOLD 

old NTM (initialized to O) 

3 

NU(20) 

bending frequency i = 1 to NCOLB (per rev) 

4 

NUNR(20) 

nonrotating bending frequency * 

i = 1 to NCOLB (rad/sec) 

bending mode displacement *1^^, i •= 1 to NBM, 
at radial station r = ' 

24 

ETA(2,10) 

^FA 

44 

ETA(2,10) 

^PB 

64 

ETA(2,10) 

^ROOT 

84 

ETA(2,10) 

1 

104 

ETA(2,10,11) 

(j-l)O.l, j “ 1 to 11 

124 

ETA(2,10,51) 

(j-l)Ar, j = 1 to MRM+1 

344 

ETA(2,10,30) 

r j = 1 to MRA 

bending mode slope 1 = 1 to NBM, at radial 

1364 


station r = ' 


ETAP(2,10) 

^FA 

1964 

ETAP(2,10) 

^FB 

1984 

ETAP(2,10) 

^ROOT 

2004 

ETAP(2,10) 

1 

2024 

ETAP(2,10,11) 

(j-l)O.l. j = 1 to 11 

2044 

ETAP(2,10,5l) 

(j-l)Ar, j = 1 to MRM+1 

2264 

ETAP(2,10,30) 

r., j = 1 to MRA 

J _i»» 

bending mode curvature i = 1 to NBM, 

3284 


at radial station r = ( 


ETAPP(2,10) 

^FA 

3884 

ETAPP(2,10) 

^PB 

3904 

ETAPP(2,10) 

^ROOT 

3924 

ETAPP(2,10) 

1 

3944 

ETAPP(2,10,ll) 

(j-l)O.l, j = 1 to 11 

3964 

ETAPP(2,10,5l) 

(j-l)Ar, j = 1 to MRM+1 

4184 

ETAPP(2,10,30) 

r j = 1 to MRA 
<3 

5204 
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MDICM 


ETAPH(2,10) 

bending mode slope at hinge, ^l^e) 

5804 

WT(ll) 

torsion frequency W,, i * 1 to NCX)LT+1, 
(per rev) 

control system frequency (per rev) 

5824 

WTO 

collective 

5835 

WTC 

cyclic 

5836 

WTR 

reactionless 

torsion mode displacement , 1 = 1 to HIM, 

at radial station r ■= 

5837 

ZETA(5,11) 

(j-l)O.l, j = 1 to 11 

5838 

ZETA(5,51) 

(j-l)Ar, j = 1 to MRM+1 

5893 

ZETA(5,30) 

r., j = 1 to MRA 

J / 

torsion mode slope 1 = 1 to NTM, at radial 

station r = 

6148 

ZETAP(5,11) 

(j-l)O.l, j = 1 to 11 

6298 

ZBTAP(5,5l) 

( j-l)Ar, j = 1 to MRM+1 

6353 

ZETAP(5,30) 

r^, j = 1 to MRA 

pltch/bendlng coupling Kp^^, 1 = 1 to MBM 

6608 

KPB(IO) 

6758 

KPG 

pitch/gimbal coupling Kp^ 

6768 

mi 

(rad) 

6769 

m2 

(rad) 

6770 

m3 


6771 

m4 


6772 

ms 

(rad) 

6773 
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INC1CM(4365) 

MB Inertia coefficients 1 
SB 2 
10 3 

IQ(lO) 4 
SQ(2,10) 

IQA(2,10) 3^ 
IQIXi(lO,10) 54 
I(iDQT(l0,10,4) 154 
IQDP(IO) 554 
IQDPT(10,4) 564 
IQIB(lO) 604 
IQIBT(10,4} 614 
SQDDP(10,5) 654 
SQDDPT(10,5,4) 704 
SQP(10,5) 904 
SQPT(10,5,4) 954 
IQO(IO) 1154 
IQ0DQ(2,10) 1164 
IQ0KiT(2,10,4) 1184 
IQODP 1264 
IQ0DPT(4) 1265 
IQOIB 1269 
IQ0DBT(4) 1270 
SQ0DDP(2,5) 1274 
SQ0DDT(2,5,4) 1284 
IFXO 1324 
IMXO 1325 
IP(5) 1326 
IPA(2,5) 1331 
IPAT( 2,5,4) 1341 
SP(2,5) 1381 
SPT(2,5,4) 1391 

ipddp(5,5) 1431 
ipddpt(5,5.4) 1456 

IPDDTT(5,5.4,4) 1556 
IPP(5.5) 1956 
spddq( 5 ,io) 1981 

SPDDCiT(5,10,4) 2031 
IP0(5) 2231 
spq( 5 ,io) 2236 

SPQT(5,10,4) 2286 
XAPQ(2,5,4,30) Xjjj at i = 1 to MRA 2486 

MQDQ( 10,10) Aerodynamic spring and damping 3686 
MQIB(lO) 3786 

mqp(io, 5) 3796 
MLQ( 10) 3846 
MDB 3856 
MP(5) 3857 
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INCICM 


QSZ 

QT 

MPIXl(5,10) 

mpib(5) 

mpdp(5,5) 

mpp(5,5) 

IQDQS(l0,10) Inertia coefficients, summed over 
IQDPS(IO) 

IQIBS(IO) 

SQDDFS(10,5) 

SQPS(10,5) 

IQODQS(2,10) 

IQODPS 

IQODBS 

SQ0DDS(2,5) 

IPAS(2,5) 

SPS(2,5) 

IPDDPS(5,5) 

SPDDQS(5.10) 

SPQS(5,10) 


3862 

3863 

3864 
3914 
3919 
3944 
3969 
4069 
4079 
4089 
4139 
4189 

4209 

4210 

4211 
4221 
4231 
4241 
4266 
4316 


(NBM=10, NTM-5, NBMT=4, MRA=30) 



WKV1CM(8165) 


CTOLD 

old 

1 

CMXOLD 


2 

GMYOLD 

old C„y 

3 

gamold(30,36) 

old (i = 1 to MRA, j = 1 to MPSI) 

4 

crcold(36) 

old max T, (j = 1 to MPSi) 

1084 

vind(3,30,36) 

V(r. , H'.) (i - 1 to MRA, j = 1 to MPSI) 

1120 

LAMBDA 

mean V 

tpp 

^GE “1“ (cos€/4z)^ (l. if OGE) 

4360 

FGE 

4361 

COSE 

cosC 

4362 

ZAGL 

^AGL 

4363 

vint(3,30,36) 

'Svint(’^i'‘^j) (i = 1 to MRA, j = 1 to MPSi) 
at other rotor 

4364 

vorh(3,36) 

Vp (j = 1 to MPSi), at other rotor huh 

7604 

LAMBDI 

mean V, a., at other rotor 
xn u 

7712 

vwb(3.36) 

( j = 1 to MPSi), at wlng/hody 

7713 

vht(3,36) 

( j “ 1 to MPSi), at horizontal tall 

-A n J 

7821 

vvt(3,36) 

Vy(Yj) ( j = 1 to MPSI), at vertical tail 

7929 

voff(3,36) 

"Vo(Vj) (j = 1 to MPSI), off rotor disk 

3037 

LAMB1]H(3) 

mean Sk y i at wing/hody 

3145 

LAMBIM(3) 

mean'Skjr, at horizontal tail 

3148 

LAMBDV(3) 

mean \y, at vertical tail 

3151 

LAMBDO(3) 

mean off rotor disk 

3154 

EINTW(3) 

®W ^ 

8157 

EINTH(3) 

8160 

EINTV(3) 

®V = ^V^V^S^F (-ks)(^R)/(S^R)ref 

8163 
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MNH1CM(462) 

ALF(10,6) 

complex (p * 1 to MHARMF), without Euler 

1 

DALF(10,6) 

suigle contributions 

complex (p = 1 to MHARMF) 

121 

DDALF(10,6) 

complex (p “ 1 to MHARMF) 

241 

PSIS(IO) 

complex (p = 1 to MHARMF) 

361 

TGOVR(lO) 

complex (^^govj.)pIj (P ” ^ ^ MHARMF) 

381 

TMAST(21) 

complex (n - 1 to MHARK) 

401 

ALF0(6) 

*%tatic 

443 

DDAL0(6) 

• 

®*^static 

4(*9 

ddalfo(6) 

“^static 

455 

PSISO 

^'^s^static 

461 

DPSISO 

static 

462 


o< 




Z. o< o< o< ) 

h X y z' 
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AES1CM(36720) 

STATE5(30,36,3) 

dcimax(30,36) 

DCEMAX(30,36) 

dcmmax(30,36) 

meff(3o,36,3) 

AEFF(30,36,3) 

DCLI3S(30,36) 

DCDES(30,36) 

dcmds(30,36) 

SAVE(30,36,19) 


integer parameter defining stall state for 
lift, drag, moment (initialized to zero) 

peak dynamic stall vortex loads (initialized 
to zero) 

'^^max 

°*^max 

'^max 

effective environment for lift, drag, moment 
Mach mimher 

angle of attack 

dynamic stall vortex load 
'^Sds 
°^ds 
®™ds 

section aerodynamic data 

(1) Up (11) 


(2) 

Up 


(12) 

°d 

(3) 



(13) 


(4) 

U 


(14) 

°dradial 

(5) 

o 

(deg) 

(15) 

F /ac 
x' m 

(6) 


(deg) 

( 16 ) 

F /ac 
r m 

(7) 

oC 

(deg) 

(17) 

F /ac 
z' m 

(8) 

M 


(18) 

M /ac 
m 

(9) 

cos^^ 

(19) 

F /ac 
T m 

(10) 

^c/v 




aerodynamic data at (r . , Hi) on disk, 
i=ltoMRA, j=l to^PST 


1 


3241 

4321 

5401 

6481 

9721 

12961 

14041 

15121 

16201 
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MNR1CM(1112) 

BETA(21,10) 

(0 

complex (l “ 1 to NBM, n = 0 to MHARM) 

1 

THETA(21,5) 

complex (i “ 1 to NTM, n = 0 to MHARM) 

421 

BETAG(21) 

complex (n “ 0 to MHARM) 

631 

PHI(10,16) 

complex (l = 1 to NAM, p » 1 to MHARMF) 

673 

PSID(10,6) 

complex (H>g ^^gg^pN 

(p = 1 to MHARMF) 

993 
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MNSCM( 12) 
QSSTAT(IO) 


1 


elastic 


PISTAT 

static 

11 


elastic 


PESTAT 

(Ve),tatl= 

12 


elastic 
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AEF1CM( 1548) 
FORCE (16,36) 

(^Ji 4 . j “ 1 to MfSI 

' j'laet rev, ** 

(dimension NBM+NTO+NGH) 

1 

FHUB(6,36) 

hub reactions (without rotor mass terms) 
F = ( ^Cjj/'^a, if2C^>^, 

577 


d2C|(^^/o-a, y2%y/v-a, -'^2C^v--a) 


T0RQUE(36) 


793 


SAVE(36,20) Integrated aerodynamic forces 

(i)-(io) 

^q^aero/®’® 

^Pj^aero/®° 

(17) Cjn^/^a 

(18) Gf /o-a 

X 

(19) Cf /o-a 

z 

(20) Cf /o-a 

r 


829 
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<iRlCM(ll39) 



QRTR(6) 

_3k T 

rotor generalized force, Q = c F 

1 

FHUBM(6) 

mean hub reaction 

# t / vJ / 

7 


F — ( ^ 2Cj^/i^ 3. t Y ZCy/ V" & » X 2 C,j^V“ 3. t 

'J(2Cni^9-a, 52C}([y/o-a, — ^2GQ/o-a) 


for trim 


CLS 



(wind 

axes) 


13 

GXS 


=x/»- 

(wind 

ajces) 


14 

GTS 






15 

GYS 


Cy/^ 




16 

GPS 


Gp/v“ 




17 

CT 

for 

inflow 




18 

GMX 


GMx 




19 

GHY 






20 

BETAO 

for 

trim 

Po 




21 

BETAC 


Pc 




22 

BETAS 


Ps 




23 

GAM(30,36) 

for 

inflow 

circulation 

n, (1 - 1 

to MBA, 0 “ 1 to MPSI) 

24 

GIRG(36) 


maximum circulation H 

(j = 1 to MPSI) 

1104 


- 43 - 



QBDCH(49) 

QWB(6) 

QHT(6) 

QVT(6) 

SAVE(3l) 


wing-body generalized forces 
horizontal tail generalized forces 
vertical tail generalized forces 
airframe aerodynamic data _ „ 


(1) 


(2) 


(3) 

(L/q)^ 

(4) 


(5) 

(My/q)ys 

(6) 


(7) 

(0/q)j^ 

(B) 

(Vq)HT 

(9) 

(D/q)vT 

(10) 

(L/q)^,^ 

(11) 

"^WB 

(12) 

^ WB 

(13) 

"*’ht 

(14) 

‘^VT 

(15) 

€ 

(16) 

<r- 

(17-19) 

V 

WB 

(20-22) 

•A 

V 

^HT 

(23-25) 

V 

VT 

(26-28) 

oS 

(29) 

4wB 

(30) 

'^HT 

(31) 

9vt 


ft or m 

1 

1 1 or m 

1 

„.2 2 

ft or m 


♦ 

deg 


ft/sec or m/sec 

I 

rad/sec 

dimensionless 

I 


1 

7 

13 

19 
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WG1CM(7998) 

RBR(3,36) 


1 

RBT(3,36) 


109 

MUTPP(3) 

/^tpp 

217 

DZT(144) 

prescribed wake, tip vortices 
D^(k), k = 1 to KRWG 

220 

DRT(144) 

D (k) , k = 1 to KRWG 

364 

K2T 


508 

DZSI(144) 

prescribed wake, sheet inside edge 
D^(k), k = 1 to KRWG 

509 

DRSl(li44) 

D (k), k = 1 to KRWG 

653 

K2SI 


797 

DZSC(144) 

prescribed wake, sheet outside edge 
D^(k), k = 1 to KRWG 

798 

DRS0(144) 

D (k), k = 1 to KRWG 

942 

K2S0 

4 

1086 

DFWG(3,2304) 

free wake, tip vortices 

D(n), n.= 1 to KRWG*MPSI 

1087 


n = (51 - l)KFVJG+k 

((k = 1 to KFWG), 1 to MPSI) 
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WKCICMC 120007) 


MR 

ML 

MI 

MW 

MH 

MV 

MO 

0 ( 3 , 20000 ) 


total number of points in flow field at which 
nonuniform induced velocity calculated for 
each azimuth (ML+MI+MW+MH+MV+MO) 

number of points on this rotor 

(MRL if INFLOW(l) = 1; zero otherwise) 

number of points on other rotor 
(MRL of other rotor if INFL0W(2) = 3i 
1 if INFLCW(2) = 2; zero otherwise) 

number of points on wing-body 

(1 if INFL0W(3) = 2; zero otherwise) 

number of points on horizontal tail 
(1 if INFL0W(4) = 2i zero otherwise) 

number of points on vertical tail 
(l if INFL0W(5) = 2; zero otherwise) 

number of points off rotor disk 
(l if INFL0W(6) = 1; zero otherwise) 

■^(n), n = 1 to HPSI*MR*MPSI 


GNW(3, 20000) 


Cj^(nj^), nj^ = 1 to MRG*(KNW+1)*MRL*MPSI 


1 

2 

3 

4 

5 

6 

7 

8 

60008 




X JA 



n = ((Si - l)*MR + k - l)*MPSI + j 

(((j = 1 to MPSI), k = 1 to MR), = 1 to MPSI) 


= (((Si - l)*MRL + k - l)*(KNW+l) + j-Sl +KNW)*MRG + i 

((((i = 1 to MRG), j = Jl-KNW to -R ), 
k = 1 to MRL), X = 1 to MPSI) 
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AEMNCM(78) 


Q(10) 

k = 1 to NBM 

1 

DQ(IO) 


11 

DDQ(lO) 

•• 

21 

p(5) 

Pj^, k = 1 to NTO (pq = p^+ p^) 

31 

DP(5) 


36 

DDP(5) 

\ 

41 

PD 

Pd 

46 

DPD 

Ph 

47 

DDPD 

Pd 

48 

PR 

Pr 

49 

DPR 

* 

Pr 

50 

DDPR 

•a 

Pr 

51 

BG 

I^G 

52 

DBG 

k 

53 

DDBG 

^G 

54 

ahub(6) 

(without Euler angle contributions to 

55 

DAHirB{6) 

•- /• • • a a «\ 

(='h \ “h 

6l 

ddahub(6) 


67 

PS 

% 

73 

DPS 

• 

74 

DDPS 


75 

TM 

^mast-bend 

76 

TG 

-^®govr 

77 

DTT 

®G " ®G ^^G 

78 


- 47 - 



LDHHCMC 2932) 

SAVEM(36,?8) motion at S>., j = 1 to MPSI- 1 

(refer to common block AEMNOM for contents) 


MB inertial coefficients for section loads 

SB 
10 

SQ(2,10) 

IQA(2,10) 

IQDQ(2,10) 

IQEB 

IQDP(2) 

SQDDP(2,5) 

SQP(2,5) 

IFXO 

IMX0(2) 

IPDDP(5) 

IPP(5) 

IPA(2) 

SPDKl(lO) 

SPQ(IO) 

SP(2) 

IPO 


2809 

2810 
2811 
2812 
2832 
2852 
2872 
2873 
2875 
2885 

2895 

2896 
2898 
2903 
29O8 
2910 
2920 
2930 
2932 
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FLMCM( 21928) 


A2(6400) 

A1(6400) 


1 

6401 

A 0(6400) 


12801 

B(2320) 


19201 

DOF1(80) 


21521 

NAMEX(80^ 


21601 

NAMEV(29) 


21681 

MX 


21710 

HXl 


21711 

MV 


21712 

MG 


21713 

DC'F1S(46) 

symmetric matrices 

21714 

NAMEXS(46) 


21760 

NAMEn/S(l6) 


21806 

MXS 


21822 

MX IS 


21823 

MVS 


21824 

MGS 


21825 

DC)F1A(43) 

antisymmetric matrices 

21826 

NAMEXA(43) 


21869 

NAMEVAU3) 


21912 

MXA 


21925 

MXIA 


21926 

MVA 


21927 

MGA 


21928 


vaxiables (80) 

^ ^^1 ^2 ^\ovri ^®govr2^ 


controls (29) 

^ ^R2 '^S 
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FLM1CM(4236) 



A2(30,30) 

^2 

1 

Al(30,30) 


901 

A0(30,30) 

*0 

1801 

AA2(30,6) 


2701 

AA1(30,6) 


2881 

aao( 30 , 6 ) 

% 

3061 

B(30,8) 

B 

3241 

BG(30,3) 


3481 

C2(6,30) 


3571 

C1(6.30) 

=1 

3751 

00(6,30) 

<=0 

3931 

ca2(6,6) 


4111 

CA1(6,6) 


4147 

CA0(6,6) 


4183 

DG(6,3) 


4219 


variables ( 30 )** Xp 
controls ( 8) ; 
gust(3)« g 
hub motion (6): 
hub forces ( 6) : F 
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FLMACM(912) 

A2(16,16) 1 

A1(16,16) 257 

A0(l6,l6) 513 

B(i 6,4) b 769 

BG(i6,3) \ 833 

BL(16,2) by^ 881 

variables (I 6 ): Xg 

controls (4) s Vg 
gust (3): g 

inflow(2): ( ) 
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FLINCM(477) 

MASSB 

1 

10 

2 

IQ(IO) 

3 

SQ(10,2) 

13 

IQA(10,2) 

33 

IQDQ(IO.IO) 

53 

IQDP(IO) 

153 

KIDB(IO) 

163 

SQDDP(10,5) 

173 

SQP(10,5) 

223 

IQ0DQ(10,2) 

273 

SQ0DDP(5,2) 

293 

IP(5) 

303 

IPA(5,2) 

308 

SP(5,2) 

318 

IPDDP(5,5) 

328 

IPP(5,5) 

353 

SPDDQ(5,10) 

378 

SPQ(5,10) 

428 
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FLAECM(646) 


MQU(IO) 

1 

MQDZ(IO) 

11 

MQZ(lO) 

MQL(IO) 

21 

31 

MQIB(lO) 

41 

MQiB(lO) 

51 

MQDQ(10,10) 

6l 

MQQ(IO.IO) 

161 

MQP(10,5) 

261 

MMU 

311 

MDZ 

312 

MZ 

313 

ML 

314 

MDB 

315 

MB 

316 

MDCi(lO) 

317 

MQ(lO) 

327 

MP(5) 

337 

TU 

342 

TDZ 

343 

TZ 

344 

TL 

345 

TIB 

346 

TB 

347 

TDQ(lO) 

34B 

TQ(lO) 

358 

TP(5) 

368 

HU 

373 

HDZ 

374 

HZ 

375 

HL 

376 

HIB 

377 

HB 

378 

HDCi(lO) 

379 

HQ(IO) 

389 

HP(5) 

399 

QU 

404 

QDZ 

405 

QZ 

406 

QL 

407 

QIB 

408 

QB 

409 

QDQ(IO) 

410 

QQ(IO) 

420 

QP(5) 

430 

RR 

435 

RU 

436 

RDZ 

437 

RZ 

438 
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FLAECM 


RL 

439 

RDB 

440 

RB 

441 

RDQ(IO) 

442 

RQ(IO) 

452 

RP(5) 

462 

MPU(5) 

467 

MPDZ(5) 

472 

MPZ(5) 

MPL(5) 

477 

482 

MPI®(5) 

487 

MPB(5) 

492 

MPDQ(5,10) 

497 

MPCi(5,10) 

547 

MPP(5,5) 

597 

mpdp(5,5) 

622 
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STDCM(882) 

DERIV(7,21) 1 
DRVR 1(7,21) (both rotors for flutter case) 148 
DRVR2(7,2l) 295 
DRVWB(7,2l) 442 
DRVHT(7,21) 589 
DRVVT(7,21) 736 


variables ( 21 ) : 

( ^F ^F % ’^F ^F ^F ^s 

®0 '^10 ®ls ^ ®lc 4 K 

“G ''g "o' 


equations (?): 

(L M N X Y Z Q) 
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S'mCH(340) 


A2FD(7,7) 

1 

A1FD(7,7) 

50 

A0FD(7,7) 

99 

BFD(7,19) 

148 

D0FFD(7) 

281 

CCNFD(16) 

288 

GUSFD(3) 

304 

D0F1FD(7) 

307 

NAMXFD(7) 

314 

NAMVFD(19) 

321 

MXFD 

340 


vaxiables ( 7 ) : 

( ^ yp 


controls ( 19 ) ! 


Ic Is 


0 Ic Is 


(e, 

^0 4c ^ 4p 


gust components in wind axes 
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TRANCM(62) 



QTRIM( 6) 

trim generalized force (total) 

1 

CQSTl 

trim —'J52CQ/^a (rotor #l) 

7 

CQST2 

trim -22G^>(r-a (rotor #2) 

8 

IB0DYI(7,7) 

inverse of body inertia 

9 

DCSAS 

SAS i. 

c 

58 

DSSAS 

SAS 

59 

mov 

transient governor 

60 

TIGOV 

transient governor ( 

61 

T2GOV 

transient governor ( A^go^)rotor#2 

62 
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2. SUBPROGRAM FUNCTION AND COMMUNICATION 


This section describes the functions of the subprograms that 
constitute the computer prograim. The communication of the subprograjns 
with each other is also described, in terms of the input and output 
variables. The description begins with the subprogram name, and its 
arguments. Next there is a statement of the principal function of the 
subprogram, and usually a general reference to a section in the analysis 
development. Then notes about the program content are given, induct ing 
references to sections in the analysis development as appropriate. 
Finally all the input and output variables of the subprogram are listed . 
The left-hand column gives the variable naune in the subprogram, and the 
right-hand column gives the label of the common block in which the 
variable is located. Some description of the variable may be given as 
well. Only the subprograms for rotor are described; the subprograms 
for rotor U7. have Identical functions and structure. 
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MAIN 


Name; MAIN 



Function: primary 

jot and analysis control 


General reference : 

section 5*3*5 


CPRTR2 

IDENT(4) 


TRIMCM 

ANTYPE(3) 

FILEID(4) 


TMDATA 

RESTRT 

JCASE 

TASK 

JOB 

RSWRT 

NCASES 

BUCDAT 

RDFILE 

START 


CASECM 
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TIMER 


Name: TIMER(N,I.T) 

Function: program timer 

N integer pareuneter controlling timing calculations 

0 initialize 

1 start timer 

2 stop timer 

3 print times 

other return present time 

I timer number 

1 case 

2 TRIM 

3 FLUT 

4 STAB 

5 TRAN 

6 STABL 

7 FLUTL 

8 WAKECl ,WAKEC2 

9 GE0MR1.GE0MR2 

10 RAMF 

11 MODE1,MOI®2 

12 M0TNR1.M0TNR2 

13 PERF 

14 LOAD 

T elapsed CPU time (sec) 

DEBUG integer parameter: print time T if GE 1 TMDATA 

ITIB 

IDB(23) 
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INPTN 


Name: INPTN 

Function: input for new job 

JCASE CASECM 

BLKDAT 

RIFILE 

DEBUGI integer parameter: debiog print control TMDATA 

OPREAD(lO) 

NROTOR 

IXX BDDATA 

lYY 

IZZ 

IXY 

IXZ 

lYZ 

ATILT 

FSCG 

BLCG 

WLCG 

WEIGHT 

FILEID(4) TMDATA 


MHARMF 
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INFTO 


Name : INPTO 

Function: input for old job 

RESTRT 

DEBUGI integer parauneter: debug print control 

NROTOR 

AHTYPE( 3 ) 

OPREAD(lO) 

DEBUG (25) 

NPRNTI 
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I 


CASECM 

TMDATA 



INPTAl 


Name: INPTAl 

Function: read airfoil table file 

DEBUG TMDATA 

TITLE(20) AITABL 

IDENT(4) 

NMAX 

NAB 

NA(20) 

A (20) 

NMB 

NM( 20) 

M(20) 

NRB 

R(ll) 

CLT(5000) 

cdt( 50 oo) 

cmt( 5000) 
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IHPTRl 


Naune t INPTRl 

Functions read rotor anraelist 

DEBUG TMDATA 

TITLE(20) RIDATA 

i 

TWISTl(5l) 
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INPTWl 


Namet INPTWl 


Function t read wake naunellst 


DEBUG 

TMDATA 

FACTWU 

• 

WIDATA 

• 

• 

KWGS0(4) 


KFWG 

• 

GIDATA 

• 

DQWG(2) 



- 65 - 



INPTB 


Name : INFIS 

Function I read body namelist 

EEBUG TMDATA 

TITLE(20) BDDATA 

ioFSYM(lO) 

LFTAM BADATA 

OPTINT 

ENGPOS ENDATA 

KEDAMP 
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INFTLl 


Name : IHPTLl 

Functions read loads namelist 

lEBUG 

MHARML 

0PN0IS(4) 

HVIB 

ZETAV(3.10,10) 


TMEATA 
LI DATA 


LA DATA 
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INPTF 


Name : INPTF 

Function: read flutter namelist for new job 

DEBUG TMDATA 

OPFIX)W FLDATA 

• 

f 

f 

NAMEXR(3) 
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INPTS 


Name : INPTS 

Function r read flight dynamics namelist for new job 

DEBUG 

NPRNTP 

• 

D0FPLT(21) 

OPTRAN 


TMDATA 

STDATA 


GCDATA 


CMAG(5) 



INPTT 


Names INPTT 

Fsinctioni read transient namelist for new job 

DEBUG TMDATA 

NPRNTT TNDATA 

OPLMDA 

OPTOAN GCDATA 

CMAG(5) 
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INPTG 


Name : INPTG 

Functions read flutter namelist for old joTs 

I]GBUG 

ANTYPE(4) 

NAMEXR(3) 


TMDATA 

FLDATA 
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INPTU 


Name ; INPTU 


Function; read flight dynamics namelist for old job 


DEBUG 

TMDATA 

0PPRNT(4) 

• 

■ 

STDATA 

• 

DCFPLT(21) 


OPTRAN 

• 

• 

GCDATA 

• 

CMAG(5) 
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INPTV 


Name: INPTV 

Function: read transient namelist for old job 


DEBUG 

NPRNTT 

NPRNTP 

NPRNTL 

NRSTRT 

TMAX 


TODATA 

TNDATA 
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FILEI 


Name: FILEI (NFILE.RDWRT) 

Function: read or write input file 

NFILE file unit number 

RDWRT integer parameter: 

TITLBD(20} 

TITLR1(20) 

TITLR2(20) 

TmCS(20) 

FILEID(4) 


all 

TMDATA 

all 

BDDATA 

all 

BADATA 

all 

ENDATA 

all 

LA DATA 

all 

GCDATA 

all 

TNDATA 

all 

STDATA 

all 

FLDATA 

all 

RIDATA 

all 

WIDATA 

all 

GIDATA 

all 

LI DATA 

all 

R2DATA 

all 

W2DATA 

all 

G2DATA 

all 

L2DATA 


0 to read file, 1 to write file 

BDDATA 

RIDATA 

R2rATA 

IMDATA 
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FILEJ 


Name: FILEJ (NFILE.RDWRT) 

Functions read or write trim data file 


NFILE file unit number 

RrWRT integer parameters 

MPSI 

LEVELl 

LEVBL2 

KNWl 

MRGl 

MRLl 

KFWGl 

KNW2 

MRG2 

MRL2 

KFWG2 

all 

all 

all 

all 

all 

all 

all 

all 

all 

all 

all 

all 

all 

all 

all 

all 

all 

all 

all 

all 

all 

all 

all 

all 

all 

all 

all 

all 


0 to read file, 1 to write file 

TMDATA 


WIDATA 


GIDATA 

W2DATA 


G2DATA 

TRIMCM 

BODYCM 

ENGNCM 

GUSTCM 

CONTGM 

CONVGM 

MNSCM 

qBDGM 

RTRICM 

RHICM 

MDICM 

INCICM 

WKVICM 

MNHIGM 

AESIGM 

MNRIGM 

AEFICM 

QRIGM 

RTR2GM 

RH2CM 

MD2GM 

ING2GM 

WKV2GM 

MNH2GM 

AES2GM 

MNR2GM 

AEF2GM 

QR2CM 
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FILER 


Name: FILER (RDWRT) 

Function: read or write restart file 


Restart file structure: 

1) case header record 

2) input, trim, airfoil data 

3) task header record — ID,NREC 
(id = 2 for flutter, 3 for flight 
dynamics, 4 for transient) 

4) task data (HFfEG records) 

5) repeat #3 and #4 as necessary 

6) end record — ID = 0, NREC = 0 


REWRT integer parameter: 0 to read file, 1 to write file 

RESTRT 

TITLCS(20) 

FILEID(4) 

NROTOR 

COIE 

IEENT(4) 

TITLR1(20) 

TITLR2(20) 

TITLBD(2Q) 

TITLA1(20) 

AF1ID(4) 

NMAXl 

CLT1(5000) 

CDT1(5000) 

CMT1(5000) 

TITLA2(20) 

AF2ID(4) 

NMAX2 

CLT2(5000) 

CDT2(5000j 

CMT2(5000) 


CASECM 

TMDATA 


TRIMCM 

HI DATA 
R2DATA 
BDDATA 

AITABL 


A2TABL 
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FILEF 


Name: 

FILEF (RDMRT) 


Function: read or write flutter restart file 


RDWRT 

integer parameter: 0 to read file, 1 to write file 


HROTOR 


TODATA 

OPFDAN 


FLDATA 

NBMl 


RTRICM 

NTOl 



HGMl 



NBM2 


RTR2CM 

N'm2 



NGM2 



all 


FLMC3M 

all 


STDCM 

all 


STMCM 

all 


MDICM 

all 


MD2CM 

all 


STDATA 

all 


GCDATA 
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FILES 


Name: FILBS(R]M?T) 

Function: read or write flight dynamics restart file 


RIWRT 

Integer parameter: 0 to read file, 1 to write file 


all 


STDCM 

all 


STMCM 

all 


STRATA 

all 


GGDATA 
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FILET 


Name: FILET(RrWRT,ENDREC) 

Function: read or write transient restart file 


RDWRT integer parameter: 0 to read file, 1 to write file 

ENDREC integer parameter: 0 if at start of transient record, 

1 if at end of record (required for file write only) 


IT 

YN(7) 

DYN(7) 

DDYN(7) 

MTRACE 
TRACE( 14377) 

LBVELl 

LEVEL2 

all 

all 

all 

all 

all 

all 


WORK 


TODATA 


TRANCM 
TNDATA 
GCEATA 
LI DATA 
L2DATA 
LADATA 
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FILEE 


Name I FILEE(KEY) 

Function: write eigenvalue file 

KEY Integer parameter defining case 

0 start file 

flutter, const, coeff, (FLUTL) 

1 complete 

2 symmetric 

3 antisymmetric 

flutter, periodic coeff. (FLUT) 

4 complete 

5 symmetric 

6 antisymmetric 
flight dynamics (STABL) 

7-18 6+IEQ (IEQ = equation type) 


TASK 

JCASE 


CASECM 

IEENT(4) 


TRIMCM 

CODE 


TMDATA 

lamda(6o) 

MX2 

\( constant coefficients) 

EIGVC 

lmdap(6o) 

lmdagp(6o) 

MX2P 

V (periodic coefficients) 

\ (periodic coefficients) 
c 

EIGVP 


- 80 - 


I 



I NIT 


Name: INIT 

Function; initialization 

NRCTOR TMDATA 
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INITA 


Name : INITA 

Function; initialize environment parameters 
General reference; section 2.5 


OPUNIT 

ALTMSL 

TEMP 

DENSEI 

OPIffiJNS 

DENSE 

ALTD 

DRATIO 

CSOUND 


TMDATA 


TRIMCM 
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INITC 


Name : INITC 

Function: initialize case parameters 

OPUNIT 

DEBUG 

MPSI 

HHARM(2) 

MHARMF(2) 

OPTRIM 

OPGOVT 

LEVEL2 

D0F(54) 

D0FT(8) 

VKTS 

VEL 

VTIP 

RPM 

COLL 

LATCYG 

LNGCYC 

PEDAL 

A PITCH 

AROLL 

ACLIMB 

AYAW 

RTURN 

NROTOR 

XTRIM 

CXTRIM 

THETFT 

PHIFT 

THETFP 

PSIFP 

THETAT 

PSIT 

D^/B0DY(6) 

DOMEGA 

DDZF 

VPIL0T(5) 

TGOVRl 

TGOVR2 

NBLDl 

A^TIPN 

RADIUS 

SIGMA 

GAMMAO 


TMDATA 


CONTCM 


RIDATA 
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INITC 


OMEGA 1 

RTRICM 

OMEGA2 

RTR2CM 

HMASS 

BODYCM 

TRATIO 

BDDATA 

CONFIG 


WEIGHT 


NBLD2 

R2DATA 

DRATIO 

TRTMGM 

DENSE 


GRAY 


CXTARG 


OFRTR2 


DPSI 


FSCALE 


RSCALE 


NSCALE 


ISCALE 


GSCALE 


SSCALE 


CSCALE 


C0SPSI(36) 


SINPSI(36) 


KEPSI(21,36) 
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INITRl 


Name: INITRl 

Function; initialize rotor parameters 

Normalization parameters: section 2.6 

Aerodynamic r, 2i.ri section 2.4.1 
Tip loss factor: section 2.4,5 

Linear twist: section 2.3*5 

Control system damping: section 5*1«3 

Gimbal/teeter spring and damping: sections 2.2.12, 2.2.13 

Lag damper; section 2.2.16 

DEBUG 

MPSI 

D0F(16) rotor degrees of freedom 

D0FT(4) 

LEVEL 

TRATIO 

DENSE 

CSOUND 

DRATIO 

QRTR(6) 

FHUBM(6) 

CLS 

CXS 

GTS 

CYS 

CPS 

CT 

CMX 

CMY 

BO 

BC 

BS 

CIRC(36) 

K2T 

K2SI 

K2S0 

GAMMAO 

SIGMA 

NBLADE 

RADIUS 

VTIPN 

TDAMPO 

TDAMPC 

TDAMPR 


TMDATA 


BDDATA 

TRIMCM 

QRICM 


WGICM 

RIDATA 
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INITRl 


NUGCO RIDATA 

NUGSO 

GDAMPC 

GDAMPS 

LDAMPC 

LDAMPM 

LDAMPR 

MRB 

MRM 

RAE(3l) 

MRA 

BTIP 

OPTIP 

TWISTApO) 

TWISTl(5l) 

RI(51) 

MRI 

INFLOW(6) 

LINTW 

TMISTL 

OMEKA RTRICM 

GLAG 

MLD 

DZLD 

CGS 

CGC 

NUGC 

NUGS 

CTO 

CTC 

CTR 

HTTP 

GAMMA 

CMEAN 

IB 

NBM 

NTM 

NGM 

NBMT 

RA(30) 

DRA(30) 

FTIP(30) 

CTOLD WKVIGM 

CMXCLT) 

CMYOLD 

vind(3,30,36) 

LAMBDA 
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INIRl 


VINT(3,30,36) WKVICM 

vorh(3,36) 

LAMBDI 

vwb(3,36) 

VHTf3,36) 

vvt(3,36) 

voff(3,36) 

LAMBIIW(3) 

LAMBDH(3) 

LAMBDV(3} 

LAMBD0(3) 

EINTW(3) 

EINTH(3) 

EINTV(3) 


STATE(30,36,3) AESICM 

DCLMAX 00 , 36 ) 

DCIMAX(30,36) 

DCMMAX(30,36) 

ALPHA ( 30 , 36 ) 

BETA(21,10) MNRICM 

THETA(21,5) 

BETAG(2l) 

PHI (10, 16 ) 

PSII)(10,6) 

QSSTAT(lO) MNSCM 

PISTAT 

PESTAT 

FORCE( 16 , 36 ) AEFIGM 

FHUB(6,36) 

T0RQUE(36) 

T750LD MDICM 

NBMOLD 

NTMOLD 

VGUST(3,30,36) GUSTGM 

VGUSTH(3) 
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IHITB 


Name: INITB 


Function; initialize airframe parameters 


Position of aircraft components: section 4.1.5 

Rotation matrix Rgp* section 4.1.2 

"Tt Rgp without rotations: sections 4.1.3t 4.1.5 

(for wind tunnel trim case) 

Control matrix section 4.1.6 

Aircraft inertia: section 4.2.4 

Airfreune elaistic modes: 

a) pitch/mast-bendlng coupling (KMST): section 4.2.3 

b) mode shape at hub (AfTODE): section 4.2.2 

c) meiss, spring, damping: section 4,2.4 

d) aerodynamic damping and control: section 4.2.7 


Initialization (for wind tunnel case) 


FV 


R = 

e 






- M* ( V X ) R^ - - M* V (tp X ) 
(Vx)Rg kp » -V 5 
G = - M* g (kp X ) 


1£BUG 

VEL 

D0F(16) 

airframe degrees of freedom 

TMDATA 

GRAY 

GAMMA 

SIGMA 

IB 

OMEGA 

NBLAQE 

RADI’JS 

reference rotor 

TRIMCM 

P21MR1 

0. 

RTRICM 

P21WR1 

(rad) 


P21MR2 

AVj, (rad) 

RTR2CM 

P21WR2 

-AM'ci (rad) 


ROTATl 

OPHVBl(3) 


RIDATA 
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INITB 


R0TAT2 

0PHVB2(3) 

R2DATA 

VGWBV(3j gust in velocity axes 

VGHTV(3j 

VGVTV(3) 

GUSTCM 

QWB(6) 

QBDGM 

QHT(6) 


QVT(6) 


amodei(6,io) 

• 

• 

BODYGM 

• 

VSIDE 


TITLE(20) 

• 

BDDATA 

K)FSYM(IO) 


DRGIW 

BADATA 
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INITE 


Name I INITE 

Function: initialize drive train parameters 

Engine inertia and control: sections 4.3*li 4,3,2 

Governor parameters (dimensionless): section 4,3,3 

Drive train spring constants: section 4,3,2 

issue 

OPENGN 

D0F(6) drive train degrees of freedom 

TRATIO 

NBLADE reference rotor 

IB 

ONEGA 

ENGPDS 

THRTLC 

lENG 

KMASTl 

KMAST2 

KICS 

KENG 

KPE 

KPl 

KP2 

TIE 

Til 

T12 

T2E 

T21 

T22 

QTHRTL 

lENGS 

KMIl 

KMI2 

KMR 

KMEl 

KME2 

KPGOVE 

KPGOVl 

KPG0V2 

TIGOVE 

TlGOVl 

T1GOV2 

T2GOVE 

T2GOV1 

T2GOV2 

Nm 


TMDATA 

BDDATA 

TRIMCM 

ENDATA 


ENGNCM 
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GHEKRl 


Name : CHEKHl 

Function: check for fatal errors 

MPSI 

LE^.^L 

NBLADE 

MHA 

RAE(3l) 

MRI 

RI(51) 

RROOT 
INFLOW (6) 

MRG 

NG(30) 

MRL 

NL(30) 

KNW 

RA(30) 

MRLO other rotor 


TMDATA 

RIDATA 

WIDATA 

RTRICM 

R2DATA 
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PRNTJ 


Name : PRNTJ 

Function; print job input data 

FILEID(4) 

all 

all 


TMDATA 

CASECM 

UNITNC* 
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PRNTG 


Name: PRNTC 

Function: print case input data 

JCASE 

JOB 

START 

FILEI:J(4) 

TrrLGS(20) 

CO I® 

ANTYPE(3) 

CPUNIT 

OPTRTH 

NRC^TC'R 

VKTS 

VEL 

RPM 

VTIP 

AlTTiSL 

TEMP 

OPGRND 

HAGL 

A FLAP 

OPENGN 

OPGC'VT 

HTIJRN 

LEVELl 

LE''EL2 

D0F(54) 

dcjftOR) 

MPSI 

MHARM 

M}t\RMF 

OPIENS 

IIENT(4) 

DENSE 

DRATIO 

CSOUND 

ALTD 

TITLE D( 20) 

WEIGHT 

FSCG 

WLGG 

BLCG 

CONFIG 

ATILT 


CASECM 

TMDATA 


TRIMCM 

BDDATA 
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PRNTC 


cws 

BODYCM 

NAM 


NDM 

ENGNCM 

TITLA1(20) 

AITABL 

AF1ID(4) 


TITU2(20) 

A2TABL 

AF2ID(4) 


TITLR1(20) 

RIDATA 

TYPEl 


RADUSl 


NBLDl 


SIGMAl 


INFLW1(6) 

0PHVB1(3) 


OPSTLl 


OFYAWl 


OPCMPl 


OPUSLl 


ROTATl 


HINGEl 


ELAGl 


EFLAPl 


GAMMAl 

RTRICM 

OMEGAl 


MTIPl / 


CMEANl ( 


IBl ^ 


NBMl 


NTMl 


NGMl 


NBMTl 


TITLR2(20) 

• 

• 

R2DATA 

■ 

EFLAP2 


GAMMA2 

• 

Rra2CM 

• 

NBMT2 
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PRNT 


Name t PRNT 

Function: print trim input data 

FILEID(4) 

MHARMF(2) 


TMDATA 
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PRNTRl 


Name* PRNTRl 

Function: print rotor input data 

NBM 

NTM 

NGM 

RA(30) 

DRA(30) 

ftipOo) 

TITLE(20) 

Tl-JISTI(51) 
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RTRICM 


RIDATA 


I 



PRNTWl 


Name: PRNTVIl 


Function? print wake input data 


MPSI 

TTIDATA 

LEVEL 


FACTOR 

WIDATA 

• 

KWGS0(4) 


KFWG 

GIDATA 

• 

DQWG(2) 
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PRNTB 


Name i PRNTB 

Function: print body input data 

NROTDR TMDATA 

TITLE(20) BDDATA 

DOFSYK(IO) 

LFTAW BA DATA 

■ 

» 

OPTINT 

ENGPOS ENDATA 

KEDAMP 
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PRNTF 


Name : PRNTF 


Fimction: print flutter input data 


IDENT(4) 

TRIMGM 

CONFIG 

BDDATA 

CPFLOW 

• 

FLDATA 

• 

GRISLD 



- 99 - 



PRNTS 


Name : PRNTS 

Function; print flight dynamics input data 

IDENT(4) TRIMGH 

NPRNTP STDATA 

« 

GUS(3) 
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PRNTJ 


Name : PRNTT 

Function: print transient input data 

IDENT{4) TRIMCM 

NPRNTT TNDATA 

OIWDA 
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PRNTG 


Name : PRNTG 

Function: print transient gust and control input data 

NROTOR 

OPTRAN 

CMAG(5) 
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I 


TMDATA 

GCDATA 



TRIM 


Name ; TRIM 
Function: trim 

General reference: sections 5«3*5i 5»3*1 


RESTRT 

RSWRT 

0PRTR2 

LEVSLl 

LEVEL2 

ITERl) 

ITERR 

ITERF 

NFRNTT 

NPRNTP 

NPRNTL 


CASECM 


'rRIHCH 


Tl'f'ATA 
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TRIMI 


LEVELl 

LEVEL2 


DEBUG 

CPTRIM 

cmiM 

CYTRIM 

BSTRIM 

BCTRIM 

OPTRIM 

MTRIM 

MTRIMD 

FACTOR 

ITERM 

ITERC 

DELTA 

EPTRIM 

OPGOVT 

CXTARG 

GRAY 

COUNTT 

CNTRLZ(ll) 

CWS 

KE(3) 

VXREKF(3) 

TCFE(ll,5) 

COUNTM 

GOUNTC 

NBLDl 

ROTATE 

NBLD2 

GAMMA 1 

OMEGAl 

IBl 

GAMMA 2 
OMEGA 2 
IB2 

VCNTRL(ll) 

THETFT 

PHIFT 

THETFP 

PSIFP 

THETAT 

DPSIF 

VPIL0T(5) 

TGOVRl 

TG0VR2 


wake analysis for rotor #1 and rotor #2: 
0 for uniform inflow, 1 for prescribed 
wake, 2 for free wake 




TMDATA 


TRIMCM 


BDDATA 

BODYCM 


CONVCM 

RIDATA 

R2DATA 

RTRIGM 

RTR2CM 


CONTCM 
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THIMI 


QRTR1(6) 

CLS 

CXS 

GTS 

CYS 

CPS 

BETAG 

BETAS 


Cp/lT- 

QRICM 

CQSl 

QRTR2(6) 

'V'" ■ 

QR2CM 

CQ32 

V-r - 



QWB(6) 

qht(6) 

QVT(6) 


QBDCM 
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TRIMP 


Name; TRIMP(LEVEL1 ,LEVEL2,ITER,ITERM) 
Function; print trim solution 


LEVELl wake analysis for rotor #1 and rotor #2; 

LEVEL2 0 for uniform inflow, 1 for inrescribed 

wake, 2 for free wake 

ITER iteration number 


ITERM 


maximum number of iterations 


CPTRIH 

CTTRIM 

CYTRTM 

BCTRIM 

BSTRIM 

OPTRTM 

HTRIM 

EITRIM 

CPGOVT 

CCLL 

LATCYC 

LNGCYG 

PEDAL 

A PITCH 

AYAW 

AROLL 

AGLIMB 

CXTARG 

GRAY 

GOUNTT 

OPRTR2 

NBLDl 

TYPEl 

NBLD2 

TYPE2 

GAMMAl 

OMEGAl 

IBl 

GAMMA 2 
OMEGA 2 
IB 2 

CWS 

KE(3) 

YXREKF(3) 


TMDATA 


TRIMCM 

RIDATA 

R2DATA 

RIRICM 

RTR2CM 

BODYCM 


- 107 - 



TRIMP 


VCNTRL(ll) CONTCM 

THETFT 

PHIFT 

THETFP 

PSIFP 

THETAT 

PSIT 

DPSIF 

VPILC’T(5) 

TGOVHl 

TGOVH2 


QRTR1(6) QPICM 


GLS 

CXS 

GTS 

GYS 

GPS 

BETAG 

BETAS 

GQSl 


= Cp/^ 


QRTH2(6) 

GQS2 


= Cp/c- 

QR2GM 

QWB(6) 

QHT(6) 

QVT(6) 


Q3DCM 
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FLUT 


Name : FLUT 

Function: flutter 

General reference: sections 5 « 3 * 5 » 5 * 3«6 


RSWRT 

CASECM 

RESTRT 


0PRTR2 

TRIMGM 

NBLADE 


OPFLOW 

FLDATA 

CPSYHM 


CPFDAN 


MPSIPG 


NINTPC 


NBLDFL 


A 2( 6^+00) 

• 

FLMGM 

MGA 


MXFD 

STMCM 
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FLUTO 


Name: FLUTM(PSl) 

Function: calculate flutter matrices 

General reference: section 6.3*1 

Inflow dynamics: sections 6.1.5* 2.4.3 

TT = 

TM = 

DLDZ = Si 

ZK = "Vt* Xv 

Drive train equations: section 6.2.3 

Construct flight dynamics matrices: section 5*3*3 also 

(only if rigid body degrees of freedom present) 
Symmetric/antisymmetric matrices: section 6.3*3 

PSI (for periodic coefficients) 

DEBUG 

OFBNGN 

0PRTR2 

DOFSYM(lO) 

TRATIO 

CONFIG 

NEM 

REULER(3.3) 

KE(3) 

RHUB1(3) 

RHUB2(3) 

AM0DE1(6,10) 

AM0DE2(6,10) 

KMSTCl(lO) 

KMSTSl(lO) 

KMSTC2(10) 

KMSTS2(10) 

MVXRE(3,3) 

TCFE(11,5) 

KIGOVE 

KIGOVl 

KIG0V2 
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IWDATA 

TRIMCM 

BDDATA 

BODYGM 


ENDATA 


I 



FLUTM 


GSE 

GSI 

QTHRTL 

lENG 

QEDAMP 

KMIl 

KMI2 

KMR 

KMEl 

KME2 

KPGOVE 

KPGC'Vl 

KPGCV2 

TIGC^’E 

TIGO' 1 

T1GOV2 

T2G0VS 

T2GCV1 

T2G0V2 

MENG22 

MENG33 

SENG22 

SENG33 

RADUSl 

NBLDl 

KFLHDl 

KHLMDl 

SIGflAl 

FXLKDl 

FYLMDl 

KINTHl 

KINTFl 

FMLMDl 

OMEGA 1 

NTMl 

NBHl 

NGMl 

MUXl 

MUYl 

MUZl 

GAMMAl 

IBl 

RGUST1(3,3) 

CHUB1(6,16) 

CBHUB1(3,3) 

GHUBT1(16,6) 


ENDATA 

ENGNCM 


RIDATA 


RTRICM 
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FLUTM 


RADUS2 

• 

• 


R2DATA 

• 

FMLHD2 



OMEGA 2 

• 


RTR2CM 

» 

CHUBT2(i6,6) 



KPBl(lO) 

KPGl 


MDICM 

KPB2(10) 

KPG2 


MD2CM 

TlCl 

TlSl 

T1C2 

T1S2 


CONTCM 

LAMBDl 

COSEl 


WKVIGM 

ZAGLl 

LAMBD2 

COSE2 

LAMBD2 


WKV2CM 

CTSl 

8 20^7- a 

QRICK 

CTS2 

y 2C^9-a 

QR2CM 

DERIV(7,21) 

DRVRl(7,2lj 

DRVWB(7,21) 

mVHT(7,2l) 

]»VVT(7,2l) 


STDCM 

A2FD(7.7) 

• 

• 


STMCK 

« 

MXFD 



CPFLOW 


FLDATA 


OPSYim 

NBLADE 

OPSAS 

KCSAS 

KSSAS 

TCSAS 

TSSAS 

0PTCRS(2) 

CPGRND 

KASGE 
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FLUTM 


D0F(80) 

CON(26) 

GUS(3) 

A 2(6400) 


FLDATA 

FLMCM 


MGA 

A2A(i6,16) 

« 

BLA(16,2) 

A2R1(30,30) FLMICM 

m 

I)GR1(6,3) 

A 2R 2(30, 30) FLM2CM 

DGR2(6,3) 
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FLUTB 


Name t FLUTB 



Fvinction: calculate flutter aircraft matrices 


General reference 

: section 6.2.2 


0PRTR2 


TRIMCM 

NEM 


BDDATA 

IB0DY(3,3) 

MSTAR 


BODYCM 

MVXRE(3,3) 

GMTRX(3,3) 

RFV(3,3) 

AMASS(lO) 

AIAMPS(IO) 

ASPRNG(IO) 

AnAMPA(lO) 

ACNTRL(4,10) 



DELTA 

OPRINT 


FLDATA 

dvbody(6) 

DDSF 

CNTRL(4) 


CONTCN 

gwb(3) 

gust in F axes 

GUSTGM 

GHT(3) 

gvt(3) 

QWB(6) 

qht(6) 

QVT(6) 


QBDG?^ 

A2(l6,l6) 

• 

• 

BL(i6,2) 


FLMAGH 

DRVWB(7,21) 

DR''HT(7,21) 

DRV\'T(7,21) 


STDCM 

LMDAW1(3) 

LMDAH1(3) 

LMDAVl(3j 

EINTWl(3) 


WKVIGM 

EINTHl(3) 

EINT\fl(3) 



LMDAW2(3) 

• 

• 

• 

EINTV2(3) 
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FLUTRl 


Name: FLUTRl (PSi) 

Function: calculate flutter rotor matrices 

General reference: sections 6.1.6, 6,4 


Azimuthal summations: 


^ fsl 

1 

at 

4- 

for periodic coefficients 

r' 

at 

> d T 

for constant coefficient 
appro xima tio n 
(section 6.1 .?) 

Reorder hub reactions: ^ 

equation multiplied by 2 to get ( - 

Inflow dynamics due 

to velocity perturbations: 

sections 6.1.4, 6,1.6 


PSI 


(periodic coefficients only) 


CPFLCW 

MPSICC 

NBLDFL 

KBM 

KTH 

NGM 

GAMMA 

NUGC 

NUGS 

GGC 

CGS 

CTO 

CTG 

CT!^ 

MUX 

MUY 

MUZ 

NBLD 

GSB(IO) 

ggt(5) 

KHLMDA 

KFLMDA 

NU(IO) 

WT(5) 

WTO 

FTC 

WTR 

KPB(IO) 

KPG 


FLDATA 


RTRIGM 


RIDATA 


MDIGM 
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FLUTRl 


LAMBDA 

WKVICM 

CTS ^ 2C^va. 

QRIGM 

TIC 

CONTCM 

TIS 


A2(30,30) 

• 

FLMICM 

DG(6,3) 


MASSB 

• 

FLINCH 

• 

SKi(5.10) 


MQU(IO) 

• 

FLAECM 


wdp(5,5) 



FLUTIl 


Name; FLUTIl (PSi) 

Function: calculate flutter inertia coefficients 

General reference; section 6.1.3 


PSI 

DEBUG TMDATA 

D0FT(4) 

GLAG RTRICM 

KBM 

Kin 

NBMT 

BETA(21,10) MNRICM 

ETAPH(2,10) MDICM 

mb INCICM 

m 

W(5,10,4) 

MASSBL FLINCM 


SPQL(5,10) 
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FLUTAl 


Name I FLUTAl (PSi) 

Function: calculate flutter aerodynamic coefficients 

General reference: section 6.1.4 

Perturbation section forces: without c/c^^ factor 

Aerodynaunic coefficients: FZO = , FXO = 

PSI 

DEBUG 

D0FT(4) 

MPSI 

DPSI 

MRA 

GH0RD(30) 

XA(30) 

XAC(30) 

CPCOMP 

OPYAW 

OPSTLL 

RFA 

RA(30) 

DRA(30) 

GMEAN 

FTIP(30) 

NBMT 

KBM 

Km 

MTIP 

MUX 

MUY 

MUZ 

ETA(2, 10 , 30 ) bending modes at r . , i = 1 to MRA 
ETAP(2,10,30) 

ETAPP(2,10,30) 

ZETA(5,30) torsion modes at r . , i = 1 to MRA 

ZETAPb.30) 

DELI 

DEL2 

DEL3 

DEL4 

DEL5 

DALPHA 

DMAGH 

OPUSLD 


TMDATA 

TRIHGM 

RIDATA 

P.TRIGM 


MDICM 

FLDATA 
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FLUTAl 


BETA(21,10) MNRICM 

DGLDS(30,36) AESICM 

DCDDS(30,36j 

DCMI]S(30,36) 

SAVE(30,36,19) 

XAPQ(2,5,^.30) INCICM 

MQU(lO) FLAECM 


™( 5 , 5 ) 
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FLUTL 


Name : FLUTL( ID , A2 , A 1 , AO ,B , MX , MX 1 , MG , DOFl , NAMEX , NAMElO 

Function: amlyze flutter constant coefficient linear equations 

Vibration point location: sections 4,l,3f 4.1.5 


ID 


A2(MX*MX} 

Al(MX*MX) 

A0(HX*MX) 

b(mx*mv) 

MX 


MXl 

MV 

MG 

DOFl(HX) 
NAMEX (MX) 

namev(mv) 

VELF(3) 

GRAV 

OMEGA 

RADIUS 

FSGG 

BLGG 

WLCG 

NSM 


THETFT 

PHIFT 

THBTAT 

PSIT 

ANTYPE(4) 

NAMEXR(3) 


problem identification: 1 for complete dynamics, 

2 for symmetric, 3 for antisymmetric 

coefficient matrices 


control matrix 

number of degrees of freedom 

number of first order degrees of freedom 

number of controls 

number of gust components 

integer vector designating first order degrees 
of freedom 

vector of variable names 
vector of control names 

BODYG'i 

TRIMCM 

reference rotor 
reference rotor 


BDDATA 


CGNTCM 


FLDATA 
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STAB 


Name : STAB 

Function: flight dynamics 

General reference: sections 5 * 3 * 5 » 5 * 3*3 

RES'fRT CASECM 

RSWRT 
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STABM 


Name j STABM 

Function; calculate flight dynamics stability derivatives and matrices 
General reference: section 5 *3 *3 

Print during stability derivative calculations; 

a) increment: 1st number dimensionless, 2nd number dimensional 

b; motion and controls: 1st number dimensionless, 2nd number 

dimensional 

1^ angular velocity = deg/sec 

2) linear velocity, gust velocity = ft/sec or m/sec 

3) Vs = rpm 

4) z-p = ft/sec^ or m/sec^ 

5) controls = deg 

c) generalized forces: moments and forces in i$2C/r'a form 

(rotor #1 parameters, body axes); 
torque in -iG-Za-a form (rotor #1 
parameters) ^ 


HPSI 

LEVELl 

LEVEL2 

DEBUG 

0PRTR2 

LSCALE 

FSCALE 

WBLDl 

MRAl 

TYPEl 

IBl 

CHUB1(6,16) 

CHUBT1(16,6) 

OMEGAl 

NBLD2 

MRA2 

TYPE2 

IB2 

GHUB2(6,16) 

CHUBT2(16,6) 

OMEGA 2 

IBODY(3,3) 

MSTAR 

MVXREp,3) 

GMTRX(3,31 

TGFE(ll,5) 
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mOATA 

TRIMCM 

RIDATA 

RTRICM 

R2DATA 

RTR2CM 

BODYCM 
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STABM 


CONFIG 

BDDATA 

QRTR1(6) 

QRICM 

CQSl 


QRTR2(6) 

QR2CM 

CQS2 -'^2C^/»-a 


101 

INCICM 

102 

INC2CM 

IRSTAR 

ENGNCM 

QTHRTL 


CiEDAMP 


KPGOVE 


KFGOVl 


KPG0V2 


KIGOVE 

ENDATA 

KIGOVl 


KIG0V2 


NFRNTP 

STDATA 

KPRNTL 


I'lERS 


OPLMHA 


DELTA 


DOF(7) 


con(i6) 


GUS(3) 


VGWBV(3) 

GUSTCM 

VGHTV(3) 


VGVTV(3) 


VGRTRlb.30,36) 


VGRTR2(3,30,36) 


VGHUB1(3) 


VGHUB2(3) 


■/GNTP.L( 11 ) 

CCNTGM 

r//B0Dy(6) 


rO^'^EGA 


DDZl- 


qwb(6) 

QBDCM 

QH?(6) 


QVT(6) 


DERIV(7,2l) 

STDCM 

• 

DRVVT(7,21) 


azfd(7.7) 

STI1CM 


M/.FD 
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STABD 


Name : STABD 

Function: print stability derivatives 

General reference: section 5*3 *3 



Dimensions i 


rotor coefficient form, M*X = y2C/«-a 
stability derivative form, X (acceleration) 
dimensionless or dimensional 


a) force or moment 


forces (FF) moments (FM) toixiue (FQ) 


M X form 
X form 


iNI^a^R 

n.% 






NlbU^ 


b) subscripts 

acceleration ( z ) 
angular velocity 
linear velocity 
controls 
gust velocity 


sA (fa) 
rz 

rjR . (FV) 

57.3 

riR (FV) 


TASK CASECM 

D0FFD(7) STMCM 

confd(i6) 

GUSFDb) 

NAMEV(19) 

ISTAR(3,3) BODYCM 

MSTAR 

IRSTAR ENGNCH 

NBLADE reference rotor TRIMCM 

IB 

OMEGA 

RADIUS 

CPPRNT(4) STDATA 

DRVR1(7,21) STDGM 

DRVR2(7,21) 

MiVWB(7,2l) 

DRVHT(7,2l) 

DRVVT(7,2l) 
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STABE 


Naime : STABE 



Function : 

calculate flight dynamics eqxiations 


DEBUG 



TMDATA 

OMEGA 


reference rotor 

TRIMCM 

EQTYEE(12) 

KCSAS 

KSSAS 

TCSAS 

TSSAS 



STDATA 

AZ^D(49) 

• 



STT'TGM 

• 

MXFD 




CPSYMM 

OPSASF 



FLDATA 

TASK 



CASECM 
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STABL 


Name: STABL(IEQ,A2,A1 .AO.B.MX ,MX1 ,MV,MG,D0F1 ,NAMEX,NAMEV,H)F,C0N) 

Function: analyze flight dynamics linear equations 

Vibration point location: sections 4.1.3i 4,1,5 

Numerical Integration of transient: sections 5*3*2, 5*3«3 

(see also program TRAN) 


lEQ 

A2(MX*MX) 

ai(mx*mx) 

A0(MX*MX) 

b(mx*mv) 

MX 

MXl 

MV 

MG 

DOFl(MX) 

NAMEX(MX) 

namev(mv) 

D0F(7) 

GGN(19) 


equation type identifier 
coefficient matrices 


control matrix 

number of degrees of freedom 

mmber of first order degrees of freedom 

number of controls 

number of gust components 

integer vector designating first order degrees of 
freedom 

vector of variable names 
vector of control names 

integer vector designating degrees of freedom used 
integer vector designating controls used 


OMEGA reference rotor 

RADIUS 

GRAY 

VELF(3) 

VGHUB1(3) 

VPTRAN(5) 

FSCG 

WLCG 

BLCG 

THETFT 

PHIFT 

THETAT 

PSIT 

DOMEGA 

NPRNTT 

D0FPLT(2l) 


TRIMCM 

BODYCM 

GUSTCM 

BDDATA 

CCNTCM 

STDATA 
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STABP 


Name t STABP( TIM , IT , YN , DYN , DDYN , DOF ) 

Function: print flight dynamics transient solution 

General reference: section 5*3*3 


Print during numerical integration (in STABL) : 

a) controls in deg 

b) gust velocity: 1st number dimensionless, 2nd ntmber 

dimensional 

c) aircraft motion: 1st number dimensionless, 2nd mmiber 

dimensional 


1) displacement = deg, ft or m 

2) velocity = deg/sec, ft/sec or m/sec 

3) acceleration = deg/sec , g 

4) inertial axes = deg/sec, g 


AANG = to == 


ALIN = a 


body 



— ( V X ) 



TIM 

time (dimensionless) 

IT 

time count 


YN(7) 

^F F F 

*F % 

DYN(7) 


h Yp h 

DDYN{7) 

( =4 ®F 


D0F(7) 

integer vector: 0 if degree of freedom not used 

GRAV 


TRIMGM 

LSCALE 



FSCALE 



TSTEP 


STDATA 

TMAX 



NPRNTT 
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STABP 


VGHUB1(3) GUSTCM 

VPrRAN(5) 

MSTAR BODYCM 

MVXRE(3,3) 

reulerO,3) 
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TRAN 


Name: TRAN 



Function: transient 


General reference: sections 5*3*5t 5*3*2 


RESTRT 

RSWRT 


CASECM 

LE1/EL1 

LET'ELZ 


TMDATA 

DVB0DY(6) 

D'^^MEGA 


CONTCM 

MVXRE(3,3) 

MSTAR 

IBODY(3,3) 


BODY CM 

ONEGA 

reference rotor 

TRIMCM 

QRTR1(6) 


QRICM 

CQSl 


QRTR2(6) 

-y2Cjj/«-a 

qR2GH 

CQS2 


QWB(6^ 

qht(6) 

QBDCM 

QVT(6) 



QTRIM(6) 

CQ3T1 

CQST2 

IB0DYI(7,7) 


TRANGM 

NPRNTT 

NPRNTP 

NPRNTL 

NRSTRT 

TMAX 

TSTEP 

OPPLOT 

D0FPLT(21) 

DOF(7) 


TNDATA 

101 


INC 1 CM 

102 


INC2CM 

CHUB1(6,16) 

GHUBT1(16,6) 

OMEGAl 

IBl 


RTRICM 
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TRAN 


CHUB2(6,16) RTR2CM 

CHUBT2(16,6) 

0MEGA2 

IB2 

nbldi ridata 

NBLD2 R2DATA 


IRSTAR 


ENGNCM 
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TRANI 


Name: TRANI(Y,DY,DDY) 


Function: 

calculate transient acceleration for 

numerical integration 

General reference: section 5*3«2 


Y(7) 

( 'I’f ^F *F ^F *F 


DY(7) 

(’+F ®F > ^F ^F V 


DDY(7) 

( ^F *^F ^F ^ ^F ®F *^s^ 


LEVBLl 

LEVEL2 

I®BUG 


TMDATA 

0PRTR2 


TRIMCM 

MVXREp,3) 

GMTRX(3,3) 


BODYCM 

TCFE(11,5) 

CNTRLZ(ll) 


BDDATA 

QTHRTL 

QEIAMP 

KPGOVE 

KPGOVl 

KPG0V2 


ENGNCM 

KIGOVE 

KIGOVl 

KIG0V2 


ENDATA 

IBl 

OMEGAl 


RTRICM 

NBLDl 


RIDATA 

IB2 

OMEGA2 


RTR2CM 

NBLD2 


R2DATA 

QBTO1(6) 

CQSl 

- ^2C|^/v-a 

QPICM 

QRTR2(6) 

CQS2 

- 52Co/o-a 

QR2CM 

QWBf6} 


QBDCM 

QHT(6j 



QVT(6) 
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TRANI 


D0F(?) TNDATA 

OPSAS 

KCSAS 

KSSAS 

TCSAS 

TSSAS 

ITERT 

OPLMm 

qtrim(6) trangm 

CQSTl 

CQST2 

IBODYlC?,?) 

DCSAS 
D6SAS 
TTGOV 
TIGOV 
T2G0V 

VPTRAN(5) 

vcmL(ii) 

DVB0DY(6) 

DOMEGA 
DDZF 

VPIL0T(5) 

TGOVRl 
TG0VR2 


GUSTGM 

OONTCM 
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TRANP 


Name; TRANP(TIM,IT,YN,DYN,DDYN) 

Function: print transient solution 

General reference; section 5*3*2 

Print notes; 

a) controls in deg 

b) gust velocity dimensionsal 

c) aircraft motion: 1st number dimensionless, 2nd number 

dimensional 

1) displacement = deg, ft or m 

2 ) velocity = deg/sec, ft/sec or m/sec 

3 ) acceleration = deg/sec , g 

4) inertial axes = deg/sec, g 

d) generalized forces; moments and forces in ll2C/^a form 

(rotor #1 pareuneters, body axes); 
torque in - yc_/cr-a form (rotor #1 
parauneters) ^ 



TIM 

time (dimensionless) 



IT 

time count 



YN(7) 

DYN(7) 

DDYN(7) 

( *^F "^F 
^F 

( ^F 

*F ^F ^ 
'♦'f ^ ^F ^ 
^F *F ^F 

'i's) 


LBVELl 

LEVEL2 




TMDATA 

FSGALE 

LSCALE 

GRAY 

0PR'ra2 




TRIMCM 
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TRANP 


ITBRT 

OFLHDA 

TSTEP 

TMAX 

MSTAR 

REULER(3,3) 

MVXREp.3) 

GMTRX(3,3) 

QTHRTL 

QEDAMP 

VGWBVf3) 

VGHTV(3) 

VGVTV(3) 

VGHUB1?3) 

VGHUB2(3) 

VPmN(5) 

NBLDl 

TYPEl 

IBI 

OMEGAl 

NBLD2 

TYPE2 

IB2 

OMEGA2 

QP'mi(6) 

CQSl 

QR'IR2(6) 

CQ32 

QVB(6) 

QHT(6) 

QVT(6) 

VCNTRL(ll) 

VPILOT(5) 

TGOVRl 

TG0VR2 

QTRIM(6) 

CQSTl 

CQ5T2 

DCSAS 

DSSAS 

TTGOV 

TIGOV 

T2GOV 


- ^ 2Cjj/i“ a 

- ^ 2Cjj/r a 


TNDATA 

BODYCM 

ENGNCN 

GUSTGM 

RIDATA 

RTRICM 

R2DATA 

RTR2CM 

QPICM 

QR2CH 

CONTCM 

TRANCM 
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TRANC 


Name: TRANC ( TIM) 


Function: 

calculate transient gust and control 


General reference: section 


TIM 

time (dimensionless) 


VELF 

MPSI 

V/flR 

TMDATA 

OMEGA 

RADIUS 

C0SPSI(36 

sinpsi(36 

0PRTR2 

reference rotor 

1 

TRIMCM 

RA1(30) 


RTRICM 

RA2(30) 


RTR2CM 

RWB(3) 

RHT(3) 

rvt(3) 

HFV(3,3) 


BODYCM 

RSF1(3,3) 



RSF2(3.3) 

RHUB1(3) 

RHUB2(3) 



MRAl 

ROTATl 


RIDATA 

MRA2 

R0TAT2 


R2DATA 

VGWBV(3) 

VGHTVp) 

VGVTV(3) 

gust in wind axes 

GUSTCM 

VGRTR1(3,30,36) 


VGRTR2(3,30,36) 

VGHUB1(3) 

VGHUB2(3) 

VPIRAN(5) 


OPTRAN 

• 

• 

• 

CMAG(5) 


GCDATA 
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CX)NTRL 


Name* OONm(T, PERIOD, C) 

Functions calculate transient control time history 
General references section 5*3»^ 


Calculates s 
T 

PERIOD 

C 


C(t) = i(l - cos2'n-t/T) 

time(sec) 
period T (sec) 
control C 
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GUSTU 


Namei GUSTU(T,PEBI0D,G) 

Function: calculate uniform gust time history 

General reference: section 

Calculates: G(t) = • 2(1 - cos2*Trt/T) 

T time (sec) 

PERIOD period T (sec) 

G gust G 
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GUSTC 


Namei GUSTC(XG,L,LO,G) 

Functiom calculate convected gust wave shape 
General reference: section 5»3»^ 


Calculates t 

G(Xg) = i(l - cos 2'Tr(Xg- L q)/L ) 

XG 

distance x (ft or m) 

L 

wavelength L (ft or m) 

LO 

starting position Lq (ft or m) 

G 

gust G 
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PERF 


Name; lERF 

Function: Performance 

General reference; section 5»2.1 


Opera tine condition; 

a; motion: 1st number dimensionless, 2nd number dimensional 


1 ) 

2 ) 

3 ) 

4) 

5 ) 

6 ) 


velocity = ft/sec or m/sec 

2 2 

dynamic pressure, q = Ib/ft or N/m 
weight, Cy/v~ = lb or N 
body motion = deg/sec, ft/sec or m/sec 
z = ft/sec^ or m/sec^ 




rpm 


b) body orientation and controls in deg 


Circulation convergence; 

a) tolerance, CG/S in form 

b) G/E = ratio'^error to tolerance (-^ 1 if converged) 

Motion convergence; 

a^ tolerance, BETA (etc) in deg 

b) BETA/E (etc) = ratio error to tolerance ( ■^ 1 if converged) 

Airframe performance; section 4,2.6 

a) aerodynamic loads: dimensional 

b) components ; 

1) angles in deg 

2) loads, q dimensional 

3 ) induced velocity, total velocity dimensionless 
Gust velocity; dimensionless 

System power; dimensional (HP); number in parentheses is percent total power 
a) climb power = V^W 

System efficiency parameters; 

a) gross weight, W = lb or N 2 

b) drag-rotor = = (Pj^ + P^)/Vj D/q-rotor = ^ * 

L/D-rotor = W/D^ ^ 

c) drag-total - - '’total/''* “A'*”*®' " ''totalA 

L/D-total - 

d) figure of merit - M - 1 - 
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PERF 


VEL 

TMDATA 

IlERM 


EIMOTN 


ITERC 


EPCIRC 


AFLAP 


0PRTR2 

TRIMCM 

G^AV 


SIGMA 


RADIUS 


OMEGA 


DENSE 


VELF(3) 

BODYCM 

VCLIMB 


VSIDE 


CVS 


HMASS 


NAM 


NDM 

ENGNCM 

NBMl 

RTRICM 

NTOl 


NGMl 


NBM2 

R1S2CM 

NTM2 


NGM2 


VGWB(3) gust in wind axes 

GUSTCM 

VGHT(3) 


VGVT(3) 
/ ^ \ 



VGHUB1(3) 

VGHUB2(3) 


VCNTRL(ll) 

CX)NTCM 

THETFT 


PHIFT 


THETFP 


PSIFP 


THETAT 


PSIT 


DVB0DY(6) 


DOMBGA 


DDZF 


save(31) 

Q3DCM 
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PERF 


LMDAW1(3) 

LMnAHl(3) 

LMDAV1(3) 

LMmW2(3j 

LMnAH2(3j 

LMDAV2(3) 

BIMS(IO) 


WKVICM 


WKV2CM 


CONVCM 


COUNTC 



PERFRl 


Namet IERFRl(P,PCPP,PI,PINT,K) ,PN) 

Function: calculate and print rotor performance 

General reference: section 5«2*1 


Operating condition: 


TPP 


1 0 fc 

0 1 Ps 

-f^c -ps 1 



o< 


HP 



© « o< 

Is TPP 



^ (^c^CP ” ^Pc ^IsW 
^Pb^Cp" ^Pb ” ^IcW 


Harmonics of glmbal motion: section 5*1*2 

Rotor forces and motion: 

shaft axes (-S), tip path plane axes (-T), wind axes (L or X) 
coefficient (Cx-), coefficient/solidity (Cx6-), dimensional (x-) 

Rotor power: LIDBAL “ '^ideal section 2.4,3) 


p total power 

PCPP climb and parasite power 


PI . induced power 

PINT interference power 

PO profile power 

PN non- ideal power 

OPUNIT 

VEL 

MPSI 

HHARM 

MHARMF 

lENSE 

NAH 

Nm 

T75 

TIC 

TIS 

FZ(30,36) F /ac 

ALPHA(30,36) 


TMDATA 


IBIMCM 

BODYCM 

ENGNCH 

OONTCM 

AESICM 
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FERFRl 


vind(3.30,36) 

LAMBDA 

. (due to other rotor) 
int ' 

WKVICM 

VIMT(3, 30 . 36 ) 
LAMBDI 

WKV2CM 

RADIUS 

SIGMA 

MBA 

TYI® 

NBLADE 

HINGE 


RIDATA 

MUX 

MUY 

MUZ 

OMEGA 

DRA(30) 

RA(30) 

ALFHP 

PSIHP 


RTRICM 

MTIP 

MAT 

NBM 

NTM 

NGM 

HUGO 

NUGS 



T750LD 

NU(20) 


MDICM 

ETA(2,10) 

bending mode at tip 


WT(ll) 

WTO 

WTC 

WTR 



FHUB(6) 

CLS 

cxs 

BETAO 

BETAC 

BETAS 

circ(36) 


QHICM 


BETA(21,10) 

THETA(21,5) 

BETAG(2l) 

PHI(10,16J 

PSID(10,6) 

qsstat(io) 

PISTAT 

PESTAT 


MNRICM 


MNSCM 



LOAD 


Name: L0AD(LEVEL1 ,LEVEL2) 

Function: loads, vibration, and noise 


Airframe vibration: section 5,2,8 

Vibration point location: sections 4,1,3, 4,1,5 

LEVELl wake analysis level for rotor #1 

LEVEL2 wake analysis level for rotor #2 

MHARMF(2) 

0PRTR2 

FSGALE 

LSCALB 

GRAY 

TRATIO 

FSCG 

WLCX; 

BLCG 

NBLDl 
ONEGA 1 
NBLD2 
OMEGA2 

MVXRE(3,3) 

MSTAR 

REULER(3,3) 

VELF(3) 

NAM 

THETAT 

PSIT 

PHI1(10,16) 

PHI2( 10,16) 

MVIB 

FSVIB(IO) 

WLVIB(IO) 

BLVIB(lO) 

ZETA(3,10,10) 


TMDATA 

TRIMCM 


BDDATA 


RIDATA 

RTRICM 

R2DATA 

RTR2GM 

BODYCM 


CONTCM 

MNRICM 

MNR2CM 

LA DATA 
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LOAIBl 


Kamei L0ADRl( LEVEL) 

Function: calculate and print rotor loads 

Print aerodynamics (function r and •V): , , ^ 

a) dimensionless quantities generally t angles in deg 

b) induced velocity in norerotatlng shaft axes 

c) interference induced velocity is that due to other rotor 

d) gust components in velocity axes 

Force/c ( dimensionless); 

L/C - iu (c/c^^)c, - 

D/0 - 

M/0 - - X/c.ean 

DR/0 » - ”radlal'%ean 

FZ/0 - OT/S - V'aean " 

'■Vo -V°mean 
"*/C -"a/Wn 

ra/= "VWn 

FRT/0 - Py Van 
Forces (dimensional) 

L - aectlon lift Wft or s/m 



D = section drag 

Ib/ft or N/m 


M — section pitch moment 

ft-lb/ft or m-N/m 


IB = section radial drag 

Ib/ft or N/m 


FZ = F = dT/dr 

Ib/ft or N/m 


z 

FX = F^ 

Ib/ft or N/m 


X 

MA = 

ft-lb/ft or m-N/m 


a 

FR = F^ 

Ib/ft or N/m 


r 

FRT = F^ 

Ib/ft or N/m 

Blade 

section power: section 5»2.1 



CP/S = d(Cp/T-)/dr 



P = section power 

HP/ft or HP/m 

XjEVEL 

level of wake analysis 




TMDATA 

OPUHIT 


MPSI 

4 It ef 



LOAERl 


dense 

DPSI 

oospsi(36) 

SIHPSI(36) 


TRIMCM 

TYPE 

RADIUS 

NBLAIE 

OPSTLL 

CHDRD(30) 

IHFL0W(6) 

HRA 


RIDATA 

OMEGA 

CMEAN 

RA(30) 

MUX 

MUY 

MUZ 

NBM 

NTO 

NGM 

PINTER(36) 

pburst(36) 


RTRIGM 

ETAT(2,10) 

ETA(2,10,30) 

bending mode at tip 

bending mode at i == 1 to MRA 

MDICM 

DBV 


WIDATA 

VGUST(3.30,36) 


GUSTGM 

GAM(30,36) 

CIRC(36) 


QRICM 

MHLOAD 

MALOAD 

MRLOAD 

RL0AD(20) 

NPOLAR 

MWKGMP 

MNOISE 

RANGE(IO) 

ELVATN(lO) 

azmuth(io) 

NPL0T(75) 


LI DATA 

SAVEM(36,?8) 


LEMNCM 

M0TI0N(78) 


AEMNCM 
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LOAIRl 


STATE(30,36,3) 

dclmax(30,36) 

k:i}hax(30,36) 

dcmmax(30,36) 

meff(30,36,3) 

aeff(30,36,3) 

dcldb(30,36) 

dcdi]b(30,36) 

dcmds(30,36) 

save(30,36,19) 

vind(3,30,36) 

LAMBDA 

vwb(3,36) 

vHr(3,36) 

vvt(3,36) 

voff(3.36) 

LAMBEWp) 

LAMBI»I(3) 

LAMBDVp) 

LAMBDO(3) 

vorh(3,36) 

vint(3,30,36) 

LAMBDI 


ABSICM 


WKVICM 


WKV2CM 
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LOAWIl 


Namet LOADHl 

Fvinctlont calculate and print hub and control loads 


Root loads: 


FHUBX = Cfj/v- 


HKUBX - CmVo- 

FHUBY “ 


MHOBZ - 

FHUBZ = Cf^A- 

•°l'cent^ 

Hub loads: 

FHUBH = Cjj/r- 

FHUBMX - CMx/'®' 


FHUBY >= Cy/«r- 

FHUBMY - 


FHUBT = Crj/r- 

FHUBQ - CqA- 

. T 

Harmonic auialysls : F = T ^ 

n j 1 

F. K 

J n 

Dimensional 

loads : 



root force i^nV(c/R) 

root moment = j?nV(c/R) 

hub force - Ni^X^^^(c/R) = J^(riR)^TrR^ 

hub moment = N^X1 L^^(c/r) = *(nR)^'»TR^ 


MHARM 

MPSI 

NBLAIDE! 

RADIUS 

TYPE 

CMEAN 

GAMMA 

OMEGA 

NBM 

NTO 

DENSE 

DPSI 

C0SPSI(36) 

SINPSI(36) 

MHARML 

NPL0T(75) 

SENDUR(12) for hub and control loads 

CMAT(12) 

EXMAT(12) 

KFATIG 
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TMDATA 

RIDATA 

RTRIGM 

TRIMCM 

LIDATA 


I 



LOADHl 


MPAER0(36) 

cmxa(36) 

CMZA(36) 

cfxa{36) 

CFZA(36) 

cfra(36) 

SAVBM(36,?8) 

MB 

SB 

10 

SQ(2,10) 

I(IA(2,10) 

IFXO 

IMXO 

IP(5) 

IPP(5.5) 

IP0(5) 

IQ0IXi(2,10) 


tV*=>a.rc 


summed over 




AEFICM 


LDMNCM 

INCICM 


SKi(5,10) 
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LOAPBl 


Names L0AI]61(R) 

Function: calculate and print blade section loads 

General reference: sections 5«2.2, 5*2. 3» 5«2.4 

Azimuth loop: PHIX = ?*'L 

PHIZ ” f 
T = © 


FXS-X = 


MXS-X = 

FXS-R = 

Cf^y^r 

MXS-Z - 

FXS-Z = 

OfjA- 

"TOR 

CENT "= 

“i'cenA 


(-■ 

B for shaft axes, P 

for principal axes' 

Harmonic analysis: 

F = ^ 2 F e- 

K„ 

Dimensional loauls: 



forces 


- ^JlV(c/R) 

moments 

“ {%/s.) 

-= ^nV(c/R) 


R radial station r/R 

MPSI 

MHARH 

D0FT(4) 

DENSE 

DPSI 

C0SPSI(36) 

SINPSI(36) 

TYPE 

NBLAI£ 

RADIUS 

HRA 

OMEGA 

CHEAN 

GAMMA 

RA(30) 

DRA(30) 

NBMT 

NBH 

NTM 
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TMDATA 

TRIMCM 

RIDATA 

RTRICM 


I 



LOADSl 


MHARML 


LIDATA 

sendur(6) 
cmat(6) 
exmat( 6) 

KFATIG 

NPL0T(75) 

for section loads 


ETA(2,10,30) 

DELI 

DBL2 

DEL3 

DEL4 

DEL5 

bending modes at i = 1 to MRA 

MDICM 

fxaero(30,36J 


AESICM 

FZAERO(30,36) 

F /ac 


maaero(30,36) 

H /ac 


fraero( 30,36) 

F^ac 


BETA(21,10) 

MNRICM 

MB 

• 

• 

• 

IPO 

SAVEM(36,?8) 


LDMNCM 
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LOADIl 


Name* LOADIl(R,Q,TR,ZR,EH? ,ER) 

Function; calculate inertia coefficients for section loads 
General reference; sections 5*2.2, 5*2.3, 5.2.4 


Blade pitch; 


section 2.3.5 

CS = cos© , SN = sine, TR =* ©(r) 


W = WP = WPP = x^ic) 

WXI = (zlC-xTc - x_k ) 

' 0 o I ' 


ZR = |-,(r), ER = ^.(r), EPR = 

^o^^trim’ ° ^*0^" *o^^trijn’ ^ 
WRXC = ( zj^ - x^K - , at r 


EPXIO(NBM) = ( 
CE(NBM) = 
CMR(MRM+l) = 


Xj) at r = e 



_* \ j * 

p - r; m d 


WFA = x^k), WPFA = (z^f - x^k)^ at r^^ 

X = 3^( J.), XR = 3^(r) 


R 

radial station ] 

Q(4) 

mean deflection 

TR 

pitch 0^ at r 

ZR(5) 

at r 

EPR(2,10) 

tk at r 

ER(2,10) 

DEBUG 

T75 

EFLAP 

ELAG 

XFA 

RFA 

ZFA 

RCPL 

NOPB 

MRM 

at r 


TMDATA 

C»NTCM 

RIDATA 
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LOADIl 


NBM 

NTW 

NGM 

NBMT 

MASS(5l) 

I‘nfflTA(5l) 

XI(51) 

TWIST(51) 

ETA(2,10,51) bending modes at r=(j-l)Ar, to MRM+1 

ETAP(2,10,51) 

ETAPP( 2,10,51) 

ZETA(5»51) torsion modes at r=(j-l)Ar, j=l to MRM+1 

ETAPH(2,10) 

EFA(2,10) bending modes at r = r„. 

EFAP(2,10) 

DELI 

IfflL2 

DEL3 

DEL4 

DEL5 

MB 


RTRICM 


MDICM 


LDMNCM 


IPO 
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LOADF 


Namet LOAIF(S,MPSI,K,SE,G,M, DAMAGE, SE5Q) 
Function* calculate fatigue damage 
General reference: section 5*2.9 


Input: 

s(mpsi) 

MPSI 

K 

SE 

M 

C 


vector of load j = 1 to MPSI; dimensional 
number of azimuthal stations; maximum 36 
parameter K in fatigue damage calculation 
endurance limit Sg (dimensional) 
material exponent 
material constant 

n 

S-N curve approximated by N = 

(S/Sj, - 1)" 


Output: 

DAMAGE damage fraction per rev (only calculated if > 0, 

C > 0, and M ^ 0) ^ 

SBQ equivalent j peak- to- peak load (only calculated if 

M ^ 0) 
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LOADM 


Name* LOADM(F,MPSI,FMEAN,FHPP) 

Function: calculate mean and half peak-to-peak 

Input: 

f(mpsi) 

MPSI 
Output: 

FHEAN mean load 

FHPP i peak-to-peak load 


load F . , j = 1 to MPSI 
number of azimuthal stations 
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GBOMPl 


LEVEL 

wake analysis: 1 for prescribed wake, 


MPSI 

2 for free wake geometry 

TMDATA 

TYIE 


RIDATA 

MWKGMP 


LI DATA 

JWKGMP(8j 

NWKGMP(4) 

KFW 


WIDATA 

KIM 

KNW 

KRW 

KRWG 

KFWG 


GIDATA 
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FOLRPP 


Name: POLRPP(A,MRA,RA,MPSI,ISUB,NPLOT,DA,NUPP) 

Function: printer-plot of polar plot 


W-NCoL* 



A 

HRA 

RA(HRA) 

MPSI 

ISUB 

NPLOT 

DA 

NUPP 


array to be plotted 

number of radial stations 

radial stations rj^, i = 1 to MRA 

number of azimuthal stations T- = j i 
j = 1 to MPSI, 360/MPSI ^ 

first dimension of array A; positive if first subscript 
is r^, negative if first subscript is 

n; data plotted every n-th step 

plot increment: last digit of integer part of 

A/DA is plotted (if multiple of NPLOT) 

unit number for printed output 
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HISTPP 


Name s HISTPP( A , MRA , RA , MPSI , ISUB , NPU)T , NAME , NUPP) 
Functions printer-plot of azimuthal time history 


d = maximum values over azimuth 

d .03 Icl 


let c = minimum, 

1 ) d >• 0, c<-. 03 d or c c 0 , 
use h = [^max(d, lc\)3 

2 ) d > 2|c \ , c > -.03 d 

use h = ■ *] 

3 ) c < - 2 |d| , d < .03 c 
use b = 

4 ) otherwise, use mean = [^ 2 ^(c + d)^ 

and b = [ max ( mean-c , d-mean ) J 


\ 4 
1 


»0 

1 

-k 

0 

w 

I 


lo4 i«7 

.. _i 

0 


U 

\ H 


lo^ 10*1 



0 

1 M 


lol 



Vo 


mean = AM = KM * 10**NM 
b = B = K * 10**N 

to convert F to K*10^ (K = 1 to 9 ) 

a) if F = 0, then F = .99 

b) N = [log 1 f\1 

if F < 1 . , then N = N - 1 

c) K = [if I /10**n] + 1 

if K = 10, then N= N+ 1 and K= 1 
if F < 0, then K= -K 

d) F = K* 10**N 
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HISTPP 


A 

MRA 

ra(mra) 

MPSI 

ISUB 

NPLOT 

NAME 

NUPP 


array to be plotted 

secondary variable; number of values (minimum l) 

secondary variable: values r. , i = 1 to MRA; 

alphanumeric labels if NPLOT ET 0; not used if 
MRA EQ 1 

number of azimuthal stations H'. = , j = 1 to MPSI, 

= 360/mpsi J 

first dimension of array A; positive if first subscript 
is r^^, negative if first subscript is 4^. 

J 

number of values of secondary variable per plot; 
minimum 1 and maximum 3» negative for alphanumeric 
labels; not used if MRA EQ 1 

name of secondary variable, 4 characters; not used 
if MRA EQ 1 

unit number for printed output 
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NOISRl 


Name: N0ISR1(RANGE,ELVATN,AZMUTH) 

Function: calculate and print far field rotational noise 

General reference: section 5»2«10 


Calculate constants: CSIR - cos^^(l - M^) 

FT = 

FD = - M^) 

FL = -N^nsinQ /4irc - M )' 

^ ® ^ ^ 

FR = -N^Hcos©^/4irc «;(! - M )‘ 

i O I- 

FB = Nricos©^/c^(l - M ) 

i S i 


FS = N^l«j/c 

Harmonic analysis of loads: 


n 
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_ -In’V, 
F . e J 
J 


RANGE 

ELVATN 

AZMUTH 

MPSI 

OFUNIT 

DPSI 

ISiNSE 

CSOUND 

cospsi(36} 

SINPSI(36) 

OMEGA 

CMEAN 

MUX 

MUY 

MUZ 

RA(30) 

IfflA(30) 

NBLAIE 

CHORD(30) 

SIGMA 

RADIUS 

MRA 

TYPE 

AXS(30) 

0PN0ISf4) 

MHARMN(3) 

MTIMEN(3) 


range s (dimensional) 
elevation ^ (deg) 
azimuth (deg) 




2 

c 


TMDATA 

TRIMCM 


RTRIGM 


RIDATA 


LI DATA 
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NCISRl 


fxa(30,36) 


AESICM 

fza(30,36) 

FVac 


fra(30,36) 



BETAC 

BETAS 


QBICM 
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BESSEL 


Name: BESSEL(NB,XB,BJ) 

Function: calculate J Bessel function 

Input: 

order of Bessel function, n 
Xg argument of Bessel function, x 

Bessel function 


Output: 

BJ 
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RAMF 


Name: RAMF(LEVEL1 ,LEVBL2,0PLMDA) 

Function; calculate rotor/airfrarae periodic motion and forces 
General reference: section 5,1.13 


Test motion convergence: section 5.1.4 

Test circulation convergence: section 5.1.12 


LEVELl 

LEVEL2 

CPLHDA 

MPSI 

MHARM( 2) 

MHARMF(2) 

ITERM 

EPMCTN 

ITBRC 

EPCIRC 

DEBUG 

MREV 

MPSIR 

0PRTR2 

NAH 

NDM 

GMEANl 

NBMl 

NTMl 

NGMl 

GMEAN2 

NBM2 

NTM2 

NGM2 

Bl(21,10) 

Tl(21,5) 

BG1(21) 

P1(10,16) 

PS1(10,6) 

B2(21,10) 

T2(21,5) 

BG2(21) 

P2(10,16) 

PS2(10,6) 


integer parameter specifying rotor #1 euid rotor #2 
wake analysis: 0 for uniform inflow, 1 or 2 for 

nonuniform inflow 

integer parameter: 0 to suppress inflow update 

TMDATA 


TRIMGM 

BODYCM 

ENGNGM 

RTRIGM 


RTR2CM 


MNRIGM 


MNR2GM 
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RAMF 


BIHS(IO) 


CONVCM 


COUHTC 

CIRC1(36) 

CTl 

CMXl 

CMYl 

CIRC2(36) 

CT2 

CMX2 

CMY2 

SIGMAl 
SIGMA 2 


QRICM 


QR2CM 


RimTA 

R2DATA 
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MODEl 


Name : H0I£1 


Function: blade modes 


T750LD 

MDICM 

NBMOLD 


NTMOLD 


I^UG 

IMDATA 

HINGE 

RIDATA 

EFMOIE 


NBM 

RTRICM 

NTM 


T75 

CONTCM 
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MOIECl 


Najne: MOIXICI 

Fiinctlon; initialize blade mode parameters 


Linearly interpolate data for bending mode calculation: section 2.3*1 

Tip mass: section 2.2.19 . xi o -s i 

Evaluate centrifugal force for bending mode calculation: section 2.3.1 

CENT = ^^i^m d^ 

Linearly interpolate data for torsion mode calculation: section 2.3.3 

Evaluate pitch inertia and control system stiffness: sections 2. 2. 9. 5* 


1 . 


3 


MRB 

MTIP 

XITIP 

EFLAP 

ELAG 

RFA 

RADIUS 

MRM 

FTO 

FTC 

FTR 

WTIN 

VTIPN 

KTOI 

KTCI 

KTRI 

MRI 



XC(5l) 

KP2(5l) 

MASS(51) 

ITHBTA(51) 

GJ(51) 

EIXX(51) 

EIZZ(5l) 

TWIST(51) 

DEBUG 


RIDATA 


TMDATA 
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MODECl 


RraiCM 

OMEGA 

EIXXB(51} 

EIZZB(51) 

MASSB(51) 

TMISTB(5l) 

CENT(5l) 

ITHETB(5l) 

GJB(5l) 

MASSI(51) 

ITHETI(51) 

xil(5lj 

XCI(51) 

TWISTI(51) 

KP2I(51) 

IPITCH 

KTO 

KTC 

KTR 
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Name: MO DEB 1 

Function: calculate blade bending modes 

General reference: section 2.3 •! 


Blade pitch: section 2.3«5 

Calculate : 

DS = 

C = 
DC = 
B = 
A = 


i _ 

s'ti 

(' 5," Cei/n*fe“T'Vi Ax- 

»C 


Normalize eigenvector solution: using Galerkin modes from last call, 

which was at r = 1 


T75 

DEBUG 


CONTCM 

TMDATA 


NOPB 

RCPL 

KFLAP 

KLAG 

EFLAP 

ELAG 

RADIUS 

RCPLS 

TSPRNG 

RFA 

RPB 

NGOLB 

MRB 

NONROT 

HINGE 

MRA 

RROOT 

MRM 

NU(20) 

NUNR(20) 

ETA(2,10,96) 

ETAP(2,10,96) 

ETAPP(2,10,96) 

ETAPH(2,10) 


RIDATA 


MDICM 
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MOEEBl 


MASSf^l) inertial and structiiral data at RIRICM 

EIXX(5l) r = e+(j-l)Ar, r = 1 to MRB+1 

EIZZ(5l) 

TWIST(5l) 

CENT(51) 

OMEGA 

MBM 

RA(30) 
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MOIEG 


Name : MDEEG(R .EFLAP ,BLAG , NCOLB , HINGE, F , IF , DDF ) 

Function I calculate Galerkin fxinctlons for bending modes 
General reference: section 2.3»1 

R radial station r/R 

EFLAP flap hinge offset e^/H 

ELAG lag hinge offset e^/R 

NGOLB number of functions 

HINGE integer parameter: 0 for hinged blade, 1 for 

cantilever blade 

F(NCOLB) Galerkin functions 

DF(NOOLB) Galerkin functions f' 

DIF ( NCOLB) Galerkin functions f* 
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MOIEAl 


Neune : MOIEA 1 


Function: calculate articulated blade flap and lag modes 

General reference: section 2.3 #2 

Calculate: F = dr, ^ ^ v|^m dr 

DEBUG 

MRB 

EFLAP 

ELAG 

KFLAP 

KLAG 

RADIUS 

MRM 

RFA 

RPB 

MRA 

RROOT 

RA(30) 

OMEGA 

NBM 

MASS(5l) section mass at r=e + (j-l)Ar, j = 1 to MRB+1 

NU(20) 

NUNR(20) 

ETA(2,10,96) 

ETAP(2,10,96) 

ETAPP(2,10,96) 

ETAPH(2,10) 


TODATA 

RIDATA 


RTRICM 


MDICM 
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MOIETl 


Name : MODETl 

Function: calculate blade torsion modes 

General reference: section 2.3*3 

Evaluate Galerkin functions at r: x = TTCr - rp^)/(l - 

Calculate. ^ ^ V;' -Lr 

B . Si WS. -u- 

<= ' 

Normalize eigenvector solution: using Galerkin functions from last 

Iteration, which was at r = 1 


IffiBUG 

MRB 

RFA 

RADIUS 

MRM 

NOOLT 

MRA 

IPITCH 

KTO 

KTC 

KTR 

OMETA 

NTM 

RA(30) 

ITHETA(51) I© at r = r„.+( j-l)Ar, j = 1 to MRB+1 

GJ(5l) GJ at r = r'J+C j-l)Ar, j = 1 to MRB+1 

WT(ll) 

WTO 

WTC 

WTR 

ZETA(5.92) 

ZETAPb,92) 


TWDATA 

RIDATA 


RTRICM 


MDICM 
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MOIEKl 


Name : HOIEKl 

Function i calculate kinematic pitch-bending coupling 
General reference* section 2.3 


lEBUG 

T75 

PHIPL 

PHIPH 

RPH 

RPB 

XPH 

KPIN 

DEL3G 

ATAHKP(lO) 

ETA(2,10) bending modes at r™ 

ETAP(2,10) 

KPB(lO) 

KPG 

NBH 


IMDATA 

OONTCM 

RIDATA 


MDICM 


RTRICM 
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MOIEDl 


Name: HODEDl 

Function: calculate blade root geometry 

General reference: section 2.2.1 


DEBUG 

T75 

CONE 

DROOP 

SWEEP 

FIROOP 

FSWBEP 

DELI 

DEL2 

DEL3 

DEL4 

DEL5 


TMDATA 

CONTCM 

RIDATA 


MDICM 
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INRTCl 


Name: INRTCl 

Function: calculate blade inertia coefficients 

General reference: section 2. 2 119 


Blade pitch: 
Calculate : 


section 2.3 «5 

CS(MRM+l) = cos©-, SN(MRM+l) = sin© 
CM(MRM-fl) = 

CMR(MRM+1) = Vf 

CMRR(MRM+l) = 

CXIM(MRM+1) = Sr 

CXIRM(MRM+l) = Sr 0 

h;m(nbm,mrm+i) = v vv. 4 ^ 

DERM(NBM,MRM+l) = Sr ^ T 5'^^ 

CEPEP(NBM,NBMT,MRM-fl) = 

X(2,NTM,NBMT,MRM+l) = 3^^ 

5,. U rl ^ 

;) X = Xj^ . for k 1 and k = 0 


a) X = 

b) XH = 


XCFA = Xp at rp^ 

XCE = at e 

XIE = Xj at e 

KP2WP = kp 


DEBUG 

T75 

MRM 

NOPB 

RCPL 

RFA 

ZFA 

XFA 

ELAG 


TMDATA 

CONTCM 

RIDATA 
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INRTCl 


RADIUS 

MBLAEB 

MRA 

EFLAP 


IB 

NBM 

NTM 

NGM 

MBMT 

RA(30) 

IPITCH 

MASS(5l) 

ITHETA(5l) 

XI(51) 

XC(51) 

KP2(5l) 

TWIST(5l) 

ETA(2,10,5l) 

BTAP(2,10,5l) 

ETAPP(2,10,5l) 

ZETA(5.51) 

ZETAP(5.51) 

EFA(2,10) 

EFAP(2,10) 

ETAPH(2,10) 

DELI 

DBL2 

DBL3 

DEL4 

DEL5 

MB 


inertial data at r= ( j-l)^ikrt j = 1 to MRM+1 


bending modes at r = (j-l)^r, 
torsion modes at r = ( j- l)Ar, 
bending modes at 


j = 1 to MRM+1 
j = 1 to MRM+1 


RIDATA 


RTRICM 


INCICM 


XAPQ(2,5.^.30) 


- 177 - 



MOEEPl 


Namei MODBPl 

Functloru print blade modes 


TYPE 

HINGE 

NCOLB 

NONRGT 

NCOLT 

RCPL 

EFLAP 

BLAG 

KFLAP 

KLAG 

RCPLS 

TSFRNG 

RADIUS 


OMEGA 

NBM 

NTW 

NGM 

NUGC 

NUGS 

KTO 

KTC 

KTR 

IB 


MB 

SB 

IP(5) 

T750LD 

NU(20) 

HUNR(20) 

ETA(2,10,11) 

ETAP(2,10,ll) 

ETAPP(2,10,11) 

WT(ll) 

WTO 


WTC 

WTR 

ZETA(5,11) 

ZETAP?5.11) 

ETAPH(2,10) 

KPB(IO) 

KPG 


bending modes at r “ ( j = 1 to 11 


torsion modes at r = (j-l).l, j = 1 to 11 


RIDATA 


RTRICM 

INCICM 

MDIGM 
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MODEPl 


DfiLl 

DEL2 

DEL3 

m4 

EEL5 


MDICM 
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BODYC 


Name: BODYC 

Function: initialize airfraune parameters at trim 

Wind tunnel trim case: section 4.1.3 

r, Rgp with rotations: sections 4.1.3, 4.1.5 

Free flight trim case: section 4.1.1 

Calculate R : section 4.2.1 

Calculate RgI*Rg, -M*(Vx)R^, G, (Vx )R^kp: section 4.2.4 

Airframe gust velocity in body axes: section 4.1.4 


THETFT 

PHIFT 

PSIFP 

THETFP 

THETA T 

PSIT 

DEBUG 

VEL 

OPTRIM 

MSTAR 

MSTARG 

ISTAR(3,3) 

RSF10(3.3) 

RSF20(3,3) 

RHUB10?3) 

RHUB20(3) 

RWB0(3} 

RHT0(3j 

RVT0(3) 

ROFFO(3) 

RSF1(3,3) 

RSF2(3,3) 

RHUB1(3} 

RHUB2(3) 

RWB(3) 

RHT(3) 

RVT(3) 

R0FF(3) 

VXRBKF(3) 

MVXRE(3,3) 

GMTRX(3,3) 

IB0DY(3.3) 

REULER(3,3) 

RFV(3,3) 


COHTGM 


THDATA 


BODYCM 
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BODYC 


rfe(3,3) 

KE(3) 

VELF(3) 

VC3LIMB 

VSIIE 

VGHBV(3) 

VGHTV(3) 

VGVTvb) 

VGWBFf3) 

VGHTF(3) 

VGVTF(3) 


BODYGM 


GUSTCM 
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ENGNC 


Najne: ENGNC 

Function: Initialize drive train parameters at trim 

Engine damping: section 4.3*1 

Drive system inertia: section 5*3 

Drive system spring, damping, mass matrices: section 5*1 >9 

Drive system static elastic matrix: section 5»1»10 

Calculate C«|) : section 5*1 *5 

Calculate Cj^ : section 5 • 1 *9 


IKBUG TMDATA 

OESNGN 

0PRTR2 TRIMCM 

NBLDl RIDATA 

NBLD2 R2DATA 

IBl RIRICM 

OMEGA 1 

GAMMAl 

CD1(2) 

CPSI1(2) 

IB2 RTR2CM 

0HEGA2 

GAMMA2 

CD2(2) 

CPSI2(2) 

101 INCICM 
QTl 

QDZl 

102 INC2CM 
QT2 

QDZ2 

CQSl QPICM 

CQS2 - ^f2C^^a 

KIGOVE ENDATA 

KIGOVl 

KIG0V2 

GSE 

GSI 

KEDAMP 
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ENGNC 


QTHRTL 

lENG 

IMIl 

KMI2 

KMR 

MKEl 

KHE2 

KPGOVB 

KPGOVl 

KPGOV2 

TIGOVE 

TlGOVl 

T1G0V2 

T2GOVE 

T2G0V1 

T2G0V2 

OEDAMP 

IRSTAR 

MENG(6,6) 

SBNG(6,6) 

DENG(6,6) 

HENG0(2,2) 


ENGNCM 
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MOTNCl 


Name ; MOTNCl 

Function: initialize rotor parameters at trim 

Calculate sections 2.4.2, 4.1.2 

Calculate R^,: section 4.1.4 

Rotor gust velocity in shaft axes: section 4.1,4 

Calculate c, c: section 4.2.2 

T 

Calculate c t section 4.2.5 
Calculate section 4.1.2 


raiBUG 

MPSI 


IMDATA 

NSCALE 

ISCALE 

FSCALE 

LSCALE 


IRIMCM 

IB 

OMEGA 

MTIP 

MUX 

MUY 

MUZ 

ALFHP 

PSIHP 

MAT 

RGUST(3,3) 

CHUB(6,16) 

CBHUB(3,3) 

CHUBT(16,6) 


RTRICM 

ROTATE 

NBLABB 

RADIUS 

MRA 


RIDATA 

NEM 


BDDATA 

DVB0DY(6) 


CONTCM 

VGUSTV(3,30,36) 

gust at rotor disk, velocity axes 

GUSTCM 

vgusts(3,30,36) 

gvist at rotor disk, shaft axes 

VGUSTH(3) 

gust at rotor hub, velocity axes 



- 184 - 



MOTNCl 


VELF(3) 

RFV(3,3) 

REULER(3,3) 

RSF(3,3) 

RHUB(3) 

AMDDE(6,10) 


BODYGM 
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BODYWl 


Name: BODYMl 


Function: 

calculate airframe transfer function matrix 


General reference: section 5*1 *8 


BEBUG 


TMDATA 

D0F(l6) 

MHARMF 

airframe degrees of freedom 


FSCALE 


TRIMCM 

NBLABB 


RIDATA 

OMEGA 

(rad) 5 0. for rotor #1 

RIRICM 

DPSI21 


CHUBT(16,6) 


AMASS(lO) 

ADAMPS{10; 

ASPRNG(10 

adampa(io 

IBODY(3,3 

MVXRE(3i3 

GMTRX(3,3! 


BODYCM 

MSTAR 

NAM 



HBODY(16,6,10) 

RHICM 
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EHGHMl 


Name : ENGNMl 

Function: calculate drive train transfer function matrix 

General reference : section 5 • 1 *9 

DEBUG 

MHARMF 

D0F(6) drive train degrees of freedom 

FSCALE 

NBLADE 

OMEGA 

DPSI21 (rad); 0. for rotor #1 

CD(2) 

MENG(6,6) 

SENG(6,6) 

DENG(6,6) 

NDM 

HENG(6,10) 


TODATA 

TRIMCM 

RIDATA 

RTRIGM 

ENGNCM 

RHICM 
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WAKEUl 


Name: WAKEUl 

Function: calculate uniform wake-induced velocity 

General reference: section 2.4,3 


Lagged thrust and moment: section 5.1.12 

Vectors for aerodynamic interference: section 4,2.6 

Interference induced velocity: section 4,2.6 

DEBUG 

OPGRND 

HAGL 

MPSI 

DPSI 

C0SPSI(36) 

SINPSI(36) 

LSCALE 

FSCALE 

MRA 

RADIUS 

ROTATE 

FACTOR 

KHLMDA 

KFLMDA 

FXLMDA 

FYLMDA 

FMLMDA 

KINTH 

KINTF 

KINTWB 

KINTHT 

KINTVT 

INFL0W(6) 

RA(30) 

OMEGA 

MUX 

MUY 

MUZ 

MRAO 

RADUSO 

OMEGAO 

RSF(3.3) 

RHUB(3) 

RWB(3) 

RHT(3) 

RVT(3) 

KB(3) 


TMDATA 


TRIMCM 


RIDATA 


RTRICM 


R2DATA 

RTR2CM 

BODYGM 
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WAKEUl 


CT C„ QPlOl 

CMY Cm„ 

CMX CmJJ 

CTOLD WKVICM 

CMXOLD 

CMYOLD 

VIND(3,30,36) 

LAHBDA 

FGE 

COSE 

ZAGL 

vint(3,30,36) 

LAMBDI 

LAMBDW(3) 

LAMBDH(3) 

LAMBDV(3) 

LAMBIX)(3) 

EINTW(3) 

einth(3) 

EIHTV(3) 
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WAKENl 


Najnei WAKENl (LEVEL) 

Fvmctionj calculate non-uniform wake induced velocity 
General reference: section 3 ^l*^ 


Calculate Nijp* section 3*1 *3 


^TS^SF 


SF^ other rotor TF 


Lagged circulation: section 5*1 *12 

Interpolate induced velocity: lineaur interpolation "between inflow 

points, constant beyond first or last 
point 

Calculate mean induced velocity: TPP normal component, aurma-weighted mean 


LEVEL 

rotor wake level: 0 for vuiiform inflow (only 

replace old circulation) 


DEBUG 

MPSI 


IMDATA 

DPSI 


TRIMCM 

MRA 

ROTATE 

INFL0W(6) 


RIDATA 

RA(30) 

DRA(30) 

DP21M 

DPSI21 

(rad); 0. for rotor #1 
(rad); rotor #2 

RTRICM 

MRAO 

ROTATO 

RA0(30) 

draoOo) 

other rotor 

R2DATA 

RTR2CM 

NG(30) 

MRG 

NL(30) 

MRL 

FACTOR 

OPVXVY 

KNW 

OrRTS 


W2DATA 

NL0(30) 

MRLO 

other rotor 

W2DATA 

RSF(3,3) 

RSFO(3,3) 

other rotor 

BODYCM 
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WAKENl 


Gam(30,36) 

QPICM 

CROC 36) 


BETAC 


BETAS 


BETA 00 other rotor 

QP2GM 

BETASO 


GAMOLD(30,36) 

WKVICM 

CR00LD(36) 



vind(3,30,36) 

LAMBDA 

vint(3,30,36) 

V0RH(3,36) 

LAMBDI 

vwb(3,36) 

vht( 3,36) 

vvt(3.36) 

voff(3,36) 

lambdm(3) 

LAMBDH(3) 

LAMBDvh) 

LAHBE0(3) 

MR WKClCM 

ML 

MI 

MW 

MH 

MV 

MO 

0 ( 3 , 20000 ) 

CHW(3i 20000) 


- 191 - 



INRTMl 


Name: INRTMl 

Function: calculate rotor transfer function matrix 

General reference: section 5*1 *6 


Aerodynamic spring and damping: section 2.2.20 


DEBUG 

D0F(15) rotor bending and 

D0FT(4) 

MPSI 

HHARM 

RA(30) 

DRA(30) 

CMEAN 

MUZ 

NUGC 

NUGS 

CGC 

CGS 

GLAG 

CTO 

CTC 

CTR 

NBM 

NTM 

NGM 

NBMT 

GAMA 

KEPSI(21,36) 

HRIR(16, 16,21) 

CT C„ 

LAMBDA ^ 

BETA(21,10} 

THETA(21,5) 

BETAG(21) 

F0RCE( 16,36) 

NBLAIE 

GSB(IO) 

GST(5) 

MRA 

CH0RD(30) 

SIGMA 

XA(30) 

XAC(30) 


TMDATA 

torsion degrees of freedom 


RIRICM 


TRIMCM 

RHICM 

QPICM 

WKVICM 

MNRICM 


AEFICM 

RIDATA 
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INTRMl 


HU(20) MDICM 

ETAPH(2,10) 

KPG 

Kra(io) 

AETA(2,10,3) bending modes at r. , i = 1 to MRA 

AZETA(5f30) torsion modes at r. , i = 1 to MRA 

WT(11) 

WTO 

WTC 

wra 

MB INCICM 

■ 

xapci(2,5,4,30) 

mqdq(io,io) 

mpp(5.5) 

IQDQS(lO.lO) 

SRIS(5,10) 
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IHRTI 


Nanet INRTI(MX,H,KEEP,LMINV,MMI»V) 

Functiom calculate inverse of transfer function matrix 


MX 

h(mx*mx) 

KSSP(HK) 

LMINV(MX+l) 

mminv(mx+i) 


dimension of H 

n 

complex matrix to be inverted 

integer vector designating degrees of freedom 
to be retained; 0 for unused degrees of freedom 

scratch vector 

scratch vector 
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MOTNHl 


Naune t MOTNHl 

Function: calculate harmonics of huh motion 

General reference: sections 


DEBUG 

MHARM 

MHARMF 

GRAY 

FSCALB 

LSCALE 

RADIUS 

ROTATE 

NBLAIB 

0PHVIB(3) 

OMEGA 

CHUB(6,16) 

CBHUB(3.3) 

CPSI(2) 

DPSI21 

KMASTC(lOj 

KMASTS(IO) 

RSF(3,3) 

KE(3) 

NAM 

NDM 

DVB0DY(6) 

DOMEGA 

QSSTAT(IO) 

PISTAT 

PHI (10, 16 } 
PSID(10,2) 
THTG(lO) 

PHIO(10,16) 

PSID0(10,2) 

THTGO(IO) 

ALF(10,6) 


(rad); 0. for rotor#l 


(due to other rotor) 

C Hs 


TMDATA 


TRIMCM 


RIDATA 


RTRIGM 


BODYCM 

ENGNCM 

OONTCM 

MNSCM 

MNRICM 

MNR2CM 

MNHICM 


DPSISO 
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MOTNRl 


Name; MOTNRl (JSTART) 

Function; calculate harmonics of rotor motion 
General reference: sections ^, 1 . 6 , 5*1 *13 


Lag damper moment: section 2.2.16 

Calculate coning and tip- path plane tilt: section 3.1.3 

Calculate hub reactions: section 5*1»7 

JSTART azimuth index 

MPSI TMDATA 

MPSIR 
DEBUG 
MHARM 
MHARMF 
D0FT(4) 

NBLADG 
GAMMA 
NBM 
NTM 
NGM 
NBMT 
GLAG 
MLD 
DZLD 
CGG 
CGS 
NUGS 
NUGC 


KPB(lO) 

MDICM 

KPG 


ETAPH(2,10) 


ETATIP(2,10) bending mode at r “ 1 


BO 

(JRICM 

BC 


BS 


BETA(21,10) 

MNRICM 

THETA(21,5) 


BETAG(2l) 


DPSI 

TRIMCM 

cospsi(36) 


SINPSI(36) 


KEPSI(21,36) 


HRTR( 16 , 16 , 21 ) 

RHICM 


RIDATA 

RTRICM 
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MOTNRl 


F0RCE(16,36) 

FHUB(6,36) 

TORQUE(3o) 

SAVE(36,20) 

Q(10) 


AEFICM 


AEMNCM 


DTT 


MB 

SB 

10 

IQ(IO) 

SQ(2,10) 

IQA(2,10) 

IQO(lO) 

IFXO 


IMXO 

IP{5) 

IPP(5,5) 

IP0(5) 

XAPQf2,5,4.30) 

MQDQ(10,10) 


INCICM 


MPP(5,5) , 

IQKidO.lO) summed over q. 

• J 

SPQ(5.10) 
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MOTNBl 


Names MOTHBl(PSl) 


Functions calculate blade and hub motion 


General references section 5*1 «5 


Rigid pitch p^s section 5*1 *3 


PSI 


q(io) 

• 

• 

AEMNCM 

• 

DTT 


MHARM 

IMDATA 

MHARFM 


NBLAIS 

RIDATA 

NBH 

RTRICH 

NTM 


NGM 


KIB(IO) 

MDICM 

KPG 


T75 

CONTCM 

TIC 


TIS 


BETA(21,10) 

MNRICM 

THBTA(21,5) 


BETAG(21) 


ALF(10,6) 

• 

MNHICM 

• 

DPSISO 
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AEROFl 


Name: AEROFl( JPSI,QT,MQ,MP,QMX,CMZ,CFX,CFZ,CFR) 

Function: calculate blade aerodynamic forces 

Calcxilate XAP * ^Aj^* section 2.2.19 

Section velocity components: section 2.4.2 

Calculate U| M, ^ : section 2.4.1 

4> In rad, In deg 

Calculate ^c/V: section 2.4.7 

Calculate cos J\- : section 2.4.6 

REVFLW = 1 If just crossed reverse flow boundary 

Tip loss correction: section 2.4.5 

Section forces and pitch moment: section 2.4.1 

FX - F/ac^. FR - F/ao^. MA - M^ao^ 

Circulation: section 2«4.9 

Unsteady lift, moment, and circulation: sections 2.4. 8, 2.4*9 

LUS = L /ac, MUS = M /ac, GUS = P /ac 

US US US 

Maxlmiun circulation outboard ^Gj^ax* section 3«1.4 


JPSI 

azimuth Index 

QT(4) 

*^jtrlm 

mq(io) 

Moj^aero/®-® 

MP(5) 

Wpjjaero/®'® 

CMX 


CMZ 


CFX 

Gf^v-a 

CFZ 

Cfz/*^a 

CFR 

Q(10) 

DQ(IO) 

DDQ(IO) 

P(5)^ 

DP(5)^ 

DBP(5) 

BG 

DBG 

DDBG 

ahub(6) 

dahub(6) 

DDAHUB(6) 



AEMNCM 
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AEROFl 


PS 

DPS 

DDPS 


AEMNCM 

DEBUG 

MPSI 


IMDATA 

DPSI 

FSCALE 

C0SPSI(36) 

sihpsi(36) 


TRIMC3M 

MRA 

CH0RD(30) 

TWIST(30) 

thetzlOo) 

XA(30) 

XAC(30) 

RGMAX 

RFA 

XFA 

OPUSLD 


RIDATA 

RA(30) 

DRA(30) 

MTIP 

OMEGA 

CMEAN 

FTIP(30) 

MUX 

MUY 

MUZ 

NBM 

NTM 

NBMT 

RGUST(3,3) 

CHUB(6,16) 


RTRIGM 

XAPQ(2,4,5^30) 


INCICM 

T75 

DVB0LY(6) 


C»NTCM 

vind(3,30,36) 


WKVICM 

vint(3,30,36) 

Interference velocity from other rotor 

WKV2CM 

GAM(30,36) 

CIRC(36) 


QRICM 

save(30,36,19) 


AESICM 

VGUST(3,30,36) 

VGUSTH(3) 

gust at rotor disk, shaft axes 
gust at rotor hub, velocity axes 

GUSTGM 
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AEROFl 


ETA(2,10,30) 

ETAP( 2 , 10 , 30 ) 

ETAPP(2,10,20) 

ZETA(5,30) 

ZETAP(5t30) 

DELI 

m2 

m3 

mA 


m5 


"bending inodes at i 

torsion modes at i 


1 to MRA 


1 to MRA 


MDICM 
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AEROSl 


Name ; AEROSl ( ALPHA , DALPHA , COSYAW , MACH, JPSI , IR ,REVFL¥ , CL , CD, CM , CDR .OPTION) 
Fxinctlom calculate Llade section aerodynamic coefficients 

Corrected Mach numbers section 2.4.5 
Stall model, delayed ocs section 2.4.7 
Yawed flow, effective c<s section 2.4.6 

Calculate 2-D airfoil characteristics at effective ©c and Ms section 2.4.7 
Section characteristics corrected for yawed flow and stall delays 


sections 2.4.6, 2.4.7 

Dynamic stall vortex loads s section 2.4.7 


ALPHA 

angle of attack cx, (deg) 


DALPHA 

<*c/V 


COSYAW 

cos-A- 


MACH 

Mach number M 


JPSI 

azimuth index j 


IR 

radial station index 1 


REVFLW 

integer parameters 1 if just crossed 
boundary 

reverse flow 

CL 



CD 

CM 

cm 

OPTION 

‘^d 

‘^dradial 

integer parameters 0 for derivatives 
in flutter analysis (no dynamic stall 
and calculated data not saved) 

of coefficients 
vortex loads. 

STATE(30,36,3) 

dclmax( 30 , 36 ) 

dcdmax( 30 , 36 ) 

DCMMAX(30,36) 

meff(30,36,3) 
aeff( 30.36,3) 
DCLDS(30,36) 
DCDDSf 30 , 36 ) 
DCMES(30,36) 


AESICM 

MRA 

MCORRL(30) 

MG0RRD(30) 

M00RRM(30) 


RIDATA 
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AEROSl 


TAUL RIDATA 

TAUD 

TAUM 

AIKLAY 

AMAXNS 

PSIDS(3) 

ALFDS(3j 

ALFRE(3) 

CLDSP 

CDDSP 

CMRSP 

OPYAW 

OPSTLL 

OPCOMP 

IffiiBUG TMDATA 

MPSI 
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AEROTl 


Name : AEROTl ( ALPHA , MA CH , RADIAL .OPTION , CL , CD , CM) 

Function: interpolate airfoil tables 

General reference: section 2,4.4 


ALPHA 

MACH 

RADIAL 

OPTION 

CL 

CD 

CM 


angle of attack o< (deg) 

Mach number M 
radial station r/R 

integer parameter: if 1 calculate c 

c , , if 3 calculate c , if 4 calculat 
Q m 

‘^’■20 


°d2D 

°m2D 


NAB 

Na(20) 

A (20) 

NMB 

NM(20) 

M(20) 

NRB 

R(il) 

CLT(5000) 

CDT(5000) 

CMT(5000) 
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^ , if 2 calculate 
e all three coefficients 


AITABL 
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BODYVl 


Name; BODYVl 

Function: calculate harmonics of airframe motion 

General reference; section 5*1 *8 


DEBUG 

MPSI 

MHARMF 

tmdata 

NBLADB 

RIDATA 

NAM 

bodycm 

hbody(i6,6,io) 

RHICM 

FHUB(6,36) 

AEFICM 

PHI (10,16) 

MNRICM 

KBPSI(21,36) 

TRIMCM 
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ENGNVl 


Name* ENGNVl 


Function: calculate harmonics of drive train motion 


General reference: section 5«1»9 


lEBUG 

™data 

MHARMF 


MPSI 


NBLAIE 

RIDATA 

NDM 

ENGNCM 

T0RQUE(36) 

AEFICM 

PSID(10,6) 

MNRICM 

heng(6,io) 

RHIGM 

KEPSI(21,36) 

TRIMCM 
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MOTNFl. 


Name; MOTNFl 

Function: calculate rotor generalized forces 

General reference: section 5*1»7 

and for trim: section 5.2.1 


DEBUG 

MPSI 

SIGMA 

GAMMA 

MUX 

MUY 

MUZ 

chubt(i6,6) 

FHUB(6,36) 

FHUBM(6) 

QRTR(6) 

CLS 

CXS 

CTS 

CYS 

CPS 

CT 

CMX 

CMY 


TMDATA 

RIIATA 

RTRICM 


AEFICM 

QRICM 
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MOTNS 


Name • MOTNS 

Function; calculate static elastic motion 
General reference; section 


MBUG 

dofa(i6) 

D0FD(6) 

airfrajne degrees of freedom 
drive train degrees of freedom 

TMDATA 

0PRTR2 


TRIMCM 

CHUBT1(16,6} 

CHUBT 2 ( 16 , 6 ) 


RTRIGM 

RTR2CM 

DDALF1(6) 

DDALF2(6) 


MNHICM 

MNH2CM 

FHUBMlf6} 

FHUBM2(6) 


QRICM 

QR2CM 

ASPRKG(IO) 

ACNTPL(4,10) 

NAM 


BODYCM 

HENG0(2,2) 

NIM 


ENGNCM 

DELF 

DELE 

DELA 

DELR 


OONTCM 

MBl 

MB2 


INCICM 

INC2CM 

QSSTAT(lO) 

PISTAT 

PESTAT 


MNSCM 



BODYF 


Name : B0DYF(LEVEL1 , LEVEL2) 

Function: calculate airframe generalized forces 

General reference: section 4,2.6 


LBVELl 

LBVEL2 

DBBUG 

MPSI 

AFLAP 

GAHMA 

SIGMA 

RADIUS 

OMEGA 

0PRTR2 

VB0DY(3) 

WB0DY(3) 

I®LF 

dele 

DELA 

DELR 

DDZF 

CANTHT 

CANTVT 

REULER(3,3) 

RWB(3) 

rht(3) 

RVT(3) 

VELF(3) 

QVfB(6) 

QHT(6) 

QVT(6) 

SAVE(31) 

VIW1(3,36) 

VIH1(3.36) 

vivi(3,36) 

LMDAWlp) 

LMDAH1(3) 

LMDAVl(3) 


wake level for rotor #1 and rotor #2: 0 for 

uniform Inflow 

TMDATA 


reference rotor TRIMC3M 


Yp 5p)_ C0NTC3M 


BDDATA 

BODYCM 


QBDCM 


WKVICM 
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BODYF 


WKV2GM 


GUSTCM 
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VIW2(3,36) 

vih2(3,36) 

viv2(3,36; 

LMnAW2(3) 

LMI)AH2(3) 

LMI1AV2(3) 

GWBf3) gust in F axes 

GHT(3) 

GVT(3) 



BODYA 


Name: B 0 DYA(VWB,VHT,VVT,WWB,AFLAP,IELF,IELE,IELA,IELR,I)AWB, 

FWB , MWB , FHT ,FVT .ANGLES ) 

Function: calculate body aerodynamic forces 

General reference: section 4.2.6 


VWB(3) 

VHT(3) 

VVT(3) 

WWB(3) 

AFLAP 

DBLF 

DELE 

BELA 

BBLR 

DAWB 

FWB(3) 

MWB(3) 

FHT(2) 

FVT(2) 

ANGLES(6) 

CANTHT 

CANTVT 

LFTAW 

OPTINT 


velocity (u, v, w) at wing-body, horizontal tail, 
and vertical tall; F axes; ft/sec or m/sec 

angular velocity (p, q, r); rad/sec 
flap angle Sp (deg) 
flaperon control (rad) 

elevator control (rad) 

aileron control ^ (rad) 

OL 

rudder control (rad) 

(rad/sec) 

(D/q, Y/q. L/q)^; ft^ or 

(M^q, My/q. ” 

(D/q, L/q)g^; ft^ or 
(D/q, L/q)y^; ft^ or 

(‘=^WB' PWB' ^VT‘ ^ 

BDDATA 

BADATA 
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WAKECl 


Name: WAKECI(LEVEL) 

Fijnction: calculate influence coefficients for nonuniform inflow 

General reference: sections 3*1»3* 3*1«^ 

Calculate h for aucisynunetric wake: section J,1.6 

Ground effect parameters: sections 2.4.3, 3.1,5 

Calculate first blade/vortex intersection age and core bursting 
age: section 3,1,? 


Make age loop: 

LANDJ = i%- 1)*MR*MPSI + j 


JTEMJ = - j 

Burst/unburst core radius; section 3«1«7 
Axisymmetrlc far wake: section 3*1 *6 

Complete C and C^^ for axisymmetric geometry: section 3*1 *6 


LEVEL 

wake analysis: 0 for uniform inflow, 1 for 

prescribed wake, 2 for free wake geometry 


NBLADE 

RADIUS 

ROTATE 

RROOT 

CHORD(30) 

MRA 

INFLOW (6) 


RIDATA 

HOTATO 

RADUSO 

other rotor 

R2DATA 

OMEGA 

CMEAN 

RA(30) 

pinter(36} 

PBURST(36) 

DPSI21 

(rad) ; - for rotor #2 

RTOIGM 

OMEGAO 

other rotor 

RTR2CM 

BETAC 

BETAS 

BETASO 

BETASO 


QRICM 

QR2CM 

MPSI 

DEBUG 

DEBUGV 

cpgrnd 

HAGL 

debug print control for VTXL and VTXS 

TMDATA 
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WAKEGl 


DPSI 

LSCALE 

FSCALE 

RWB(3) 

RHT(3) 

RVT(3) 

RHUB(3) 

RHUB0(3) other rotor 

ROFF(3) 

RSF(3,3) 

RSFO(3,3) other rotor 

KE(3) 

RFE(3,3) 

K2T 

MUTPP(3) 

KNW 

KRW 

KFW 

KDW 

RRU 

FRU 

PRU 

FNW 

DVS 

DLS 

CORB(5) 

OPCORE(2) 

WKMOnL(l3) 

OPHMS(2) 

LHU 

OPHW 

OPRTS 

VELB 

DPHIB 

rffiV 

QCEBUG 

MRG 

ng( 30 ) 

MRL 

NL(30) 

MRLO other rotor 


TRIMCM 


BODYCM 


WGICM 

WILATA 


W2DATA 
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WAKECl 


MR 

ML 

MI 

MW 

MH 

MV 

MO 

0 ( 3 , 20000 ) 

CNW(3,20000) 


WKClCM 
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WAKEBl 


Name: WAKEBl(PSI .OPTION, RBR.RBT.RB) 
Function: calculate blade position 
General reference: section 3 •1*3 


PSI 

OPTION 


RBR(3) 

RBT(3) 

RB(3,30) 

MPSI 

MHARMF 

MHARM 

RFA 

ZFA 

XFA 

NBLAIX) 

RROCT 

NBM 

RA(30) 

OFWKBP(3) 

MRG 

NG(30) 

MRL 

NL(30) 

BETA(21,10) 

BETAG(2l) 

PSIS(IO) 

PSISO 

ETA(2,10,30) 

ETAR(2,10) 

ETAT(2,10) 

DELI 

DEL2 

DEL3 


(rad) 

integer parameter controlling calculation of 

if 1, at and 1; if 2, at circulation stations; 

if 3, at inflow stations 

^b ’^ROOT 
r, at tip (r = l) 

D 

IV at inflow or circulation stations 
b 

TMDATA 

RIDATA 


RTRICM 

WIDATA 


MNRICM 

MNHICM 


bending modes at r^, i = 1 to MRA 
bending modes at I’pqqt 
bending modes at tip Cr = l) 


MDICM 
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VTXL 


Name : VTXL(R1 ,R2,RP,M0IEL,0PCX)RE , CX3RE, DLS , CHORD, PSI ,OPGRND, ZAGL.RTE , 

VI, v2 , debug) 

Function; calculate vortex line segment induced velocity 
General reference; section 3*1 •? 

Calculate; SI = sVs, S2 = s^Si RMSQ = r^ 

1 c m 

Lifting surface correction; 

ANGLS = J'L(deg) 

HLS = h (-1.0 for no correction) 

RSINL = rsin-A_ , COSL = cos_A_, SINL = sinJV> 

LLL = LLS = FACTLS = 

Image element in ground effect; section 3*1 »5 

Rl(3) 

R2(3) 

RP(3) 

MODEL 

0 POORE 

CORE 
DLS 
PSI 
CHORD 
OFGRND 
ZAGL 

RTE(3,3) 

DEBUG 

Vl(3) 

V2(3) 
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(at «j>) 

r^ (at 4 *+aH^) 
rp (atP) 

integer parameter; 1 for stepped vorticity distribution, 
2 for linear vorticity distribution 

integer parameter defining vortex core type; 0 for 
distributed, 1 for concentrated vorticity 

vortex core radius r 

c 

d^g for lifting surface correction, LT 0. to suppress 

^ ; required for d^^^ ^ 0 only 

chord c at P; required for d^^ '5^. 0 only 

integer parameter; 0 for out of ground effect 

^AGL’ ground effect only 

R^g; required in ground effect only 

integer parameter; debug print if GE 3 

Av due to (at t^) 

Av due to Pp (at c^+ A.'p) 



VTXS 


Name ; VTXS (R1,R2,R3,R4,RP,M0DELT,M0 DELS , OPCORE , CCRBT , CORES , DV'S , 

OPGRND, ZAGL.RTE.MDLT.MDLS.VTl ,VT2,VS1 ,VS3 .DEBUG) 

Function: calculate vortex sheet segment induced velocity 

General reference: section 3»1»8 

Image element in ground effect: section 3 •1*5 


Rl(3) 

R2(3) 

R3(3) 

R4(3) 

RP(3) 

MODEL? 

MODELS 

OPCORE 

CORE? 

CORES 

DVS 

OPGRND 

ZAGL 

RTB(3,3) 

DEBUG 

MDLT 

MDLS 

VT1(3) 

VT2(3) 

VS1(3) 

VS3(3) 


^2 

^4 

’^P 

integer parameters defining trailed and shed vorticity 
model: 0 to omit, 1 for stepped line, 2 for linear 

line, 3 for sheet 

integer parameter defining vortex core type: 0 for 

distributed, 1 for concentrated vorticity 


for trailed vorticity (LT 0. for s/z) 

r for shed vorticity (LT 0. for t/2) 
c 

d for sheet edge test; LT 0. to suppress 

VS 

integer parameter: 0 for out of ground effect 

z. -- ; required in ground effect only 
AGlf 

R— ; required in ground effect only 

To 

integer parauneter: debug print if GE 3 

integer parameters specifying trailed and shed vorticity 
model used 


(at outside edge) 

T 2 (at outside edge) 

Av due to r. (at outside edge) 

S X 

due to inside edge) 


Av^ due to 
due to 


Av 



GEOHEl 


Name: GB0ME1(K,L, LEVEL, RWT,RWSO,RWSl) 

Function: evaluate wake geometry 

General reference: section 3 * 1 *3 


K 

k ( *= k ) 


L 



XiETV^Xj 

wake analysis: 1 for prescribed wake geometry, 

free wake geometry 

2 for 

RWT(3) 

"r at tip vortex 


RWS0(3) 

at sheet inside edge 


R¥SI(3) 

r^ at sheet outside edge 


MPSI 


TMDATA 

DPSI 


miMCM 

KRWG 


WIDATA 

KFWG 


GIDATA 

rbr(3,36} 

rbt(3,36) 

MUTPP(3) 

DZT(144) 

I»T(144) 

K2T 


WGICM 

DZSI(144) 

DRSI(144) 

K2SI 

DZS0fl44) 

DRS0(144) 

K2S0 

DFWG(3,2304) 
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GBOMRl 


Name: GEOMRl( LEVEL) 

Function: calculate wake geometry distortion 

General reference: section 3*1 *3 


Prescribed wake geometry: CTG ■= C,j,, CTOS - TW = (deg) 


LEVEL 


DEBUG 
MPSI 
DPS I 

NBLADE 

SIGMA 

TMIST(30) 

KHLMDA 

RROOT 

MRA 

LAMBDA 

LAMBDI 

KBVG 

OPRWG 

FWGT(2) 

FWGSI(2) 

FWGSO(2) 

KWGT(4) 

KWGSI(4} 

KWGS0(4) 

CT 

CIRC(36) 

BETAC 

BETAS 

HA(30) 

MUX 

MUY 

MUZ 

RBR(3,36) 


wake analysis: 1 for prescribed wake geometry, 2 for 

free wake geometry 

IMDATA 

TRIMCM 

RIDATA 


at r^, i = 1 to MRA 


interference velocity, due to other rotor 


WKVICM 

WKV2CM 

WIDATA 


QPICM 


RTRICM 


WGICM 


K2S0 
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GEOMFl 


Name : GEOMFl 

Function: calculate free waJce geometry distortion 

General reference: section 3.2 


Subprograms required: WGAM, DCALC, NWGAL, WQCAL, VSCAL, Q3VL, QCVL, QVS 


DEBUG 

integer parameter controlling debug 
print: GE 1, print D at ^ = 2 it/N each 

TMDATA 


iteration; GE 2, allow printing; GE 3. 
controlled by IPWGDB and QNGIB 


MPSI 

(maximum 24, multiple NBLAEE) 


SIGMA 

NBLAIE 


RIDATA 

phibwg(36) 

core burst age (rad) 

RTRICM 

DBV 

WIDATA 

MUTPP(3) 

DFWG(3,2304) 


WGICM 

LAMBDA 

FACTGE 


WKVICM 

LAMBDI 

interference velocity, due to other rotor 

WKV2CM 

CONING 

& (rad) 

QPICM 

CIRC(36) 

P/rA 


KFWG 

CPFWG 

IlERWG 

FACIWG 


GIDATA 

WGM0DL(2) 

RTWG(2) 



G0REHG(4) 

MRVBWG 

LIMWG 



NDMWG(36) 

IPWGIB(2) 



QWGIB 

DQWG(2) 



mi 


MDICM 

DEL2 
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MINV 


Name: MINV(A,N,D,L,M) 

Function: calculate inverse of matrix 

Input; 

a(n*n) 

N 

L(N+l) 

M(Ntl) 

Output: 

A(N*N) a - inverse 

D determinant of A; 0. if A is singular 


matrix (destroyed) 
dimension 
scratch vector 
scratch vector 


- 221 - 



MINVC 


Name: MINVC(A,N,D,L,M) 

Function: 

calculate inverse of complex matrix 

Input: 

a(n*n) 

complex matrix 

N 

dimension 

L(N+l) 

scratch vector 

M(N+1) 

scratch vector 

Output: 

a(n*n) 

complex A - inverse 

D 

complex determinant of A; 0. if A is singular 
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EIGENJ 


Name i EIGENJ ( N , NM , A , T , EVR , EVI , VEGR , VECI , INDIC , NEI ) 
Function* calculate eigenvalues and eigenvectors of matrix 

Subprograms required: SCALEM, HE!SQNt REALVE, COMPVE 


Input* 

A(N*N) 

N 

NM 

NEI 

T 


matrix A (destroyed) 
order of matrix 

actual first dimension of arrays; maximum 100 
0 to calculate only eigenvalues 
dummy argument (set to 24. in EIGENJ) 


Output* 

evr(n) 

evi(n) 

VECR(N*N) 

VECI(N*N) 

INDIC(N) 


real part of eigenvalues of A 

imaginary part of eigenvalues of A 

real part of eigenvectors of A 

imaginary part of eigenvectors of A 

if 2, no error; if 1, eigenvector not found; if 0, 
neither eigenvector nor eigenvalue found 
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lERED 


Name» EERBD(NX,NV,D0F,C0N,A2,A1,A0,B,I]0F1,I]0F0,NAMEX,NAMEV) 

Function* eliminate equations and variables from system of differential 
equations 


Input* 

NX dimension of matrices 


HV 

dof(nx) 

oon(nv) 

A2(NX*NX) 

A1(NX*NX) 

AO(NX*NX) 

B(NX*NV) 

DOFO(NX) 

DOFl(NX) 

NAMEX(NX) 

NAMEV(NV) 


dimension of matrices 

integer vector designating degrees of freedom to be 
eliminated* DOF = 0 if variable not used 

integer vector designating controls to be eliminated: 
CON = 0 if variable not used 

coefficient matrices 


control matrix 
integer vector 
integer vector 
vector of variable names 
vector of control names 


Output* 

A2 reconstructed matrices and vectors 

A1 
AO 
B 

DOFO 

LOFl 

NAMEX 

NAMEV 
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Q3TRAN 

Name* QSTRAN(MX,MX0,MX1,MV,A2,A1,A0,B0,D0F1,D0F0,NAMEX) 

Fxmctlon* quasistatic reduction of system of linear differential 
equations 

General reference: section 6.3 .2 


Input* 

A2(MX*MX} 

Al(MX*MXj 

A0(MX*MX) 

BO(HX*MV) 

EOFl(MX) 

dofo(mx) 

MX 

MXO 

MXl 

MV 

NAMEX(MX) 


coefficient matrices 


control matrix 

integer vector designating first order degrees 
of freedom: DOFl(l) = 0 for first order 

integer vector designating quasistatic variables* 
D0F0(l) = 0 for quasistatic 

number of degrees of freedom, maximum 60 

nmber of quasistatic degrees of freedom 

number of first order degrees of freedom 

number of controls , maximum 60 

vector of variables names 


Output * 

A2 reconstructed matrices and vectors 

A1 

AO 

BO 

DOFl 

NAMEX 

MX number of remaining degrees of freedom (MX-MXO) 

MXl number of remaining first order degrees of freedom 
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CSYSAN 


Namei CSYSAN(N,MX,MX1,MV,A2,A1,A0,B0,MFRBQ,FREQ,NSTEP,II)F1,FSCALE, 

nahex.namev.mfout) 

Function: analyze system of constant coefficient linear differential 

equations 

General reference: sections 7.2, 7.2,1 


N 


A2(MX*MX} 

A1(MX^»MX) 

AO(MX^m; 

BO(MX^W) 

MX 

MXl 

MV 

DOFl(MX) 

FSCALE 

namex(mx) 

nambv(mv) 

NSTEP 

NFREiQ 

freq(nfreq) 

NFOUT 


calculation control 


o 

il 

1 

2 

10 

11 

12 

eigenvalues x 

X 

X 

X 

X 

X 

eigenvectors 

X 

X 


X 

X 

check sums 


X 



X 

zeros 



X 

X 

X 


coefficient matrices 


control matrix 

number of degrees of freedom 

nmber of first order degrees of freedom 

number of controls 

(maximum MX2= 2*MX- MX1= 60; maximum MV = 60) 

Integer vector designating first order degrees of 
freedom (zero columns in AO); DOFl(l) = 0 for x^^ 
first order 

frequency scale factor IT. (in rad/sec to obtain 
frequencies in Hz and times in sec); there is no 
print of dimensional eigenvalues if FSCALE ^ 0 . 

vector of variables names 

vector of control names 

static response calculated if NS’IEP 4 0 

number of frequencies for which frequency response 
calculated; none if NFREQ ^ 0 

vector of frequencies (dimensionless) for calculation 
of frequency response 

unit number for printed output 
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CSYSAN 


Outputt 

LAMDA(MX2) 

MX2 


eigenvalues 

number of eigenvalues 

available in following common block: 

COMMON /EIGVC/LAMI]A(60),MX2 
COMPLEX LAMDA 
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EETRAN 


Name: EETRAN(A,MX,MX1,MV,A2,A1,A0,B0,D0F1,NAMEX,NAME,NP0UT) 

Function: transform equations to state variable form 

General reference: section 7*1 


Input: 

A2(MX*MX} 

A1(MX*MX) 

ao(mx*mx) 

B0(MX*MV) 

MX 

MXl 

MV 

IX)Fl(MX) 

namex(mx) 

NFOUT 


coefficient matrices 
control matrix 

number of degrees of freedom, maximum 60 
number of first order degrees of freedom 
number of controls, maximum 60 

integer vector designating first order degrees of 
freedom; IDFl(l) = 0 for Xj^ first order 

vector of variable names 

unit number for printed output 


Output: 

A(MX2*MX2) coefficient matrix 

B0(KX*MV) control matrix 

NAME(MX2) vector of variable names 

(MX2 = 2*MX - MXl) 
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SINE 


Name » SINE( W , A , ASQ ,B0 ,MX , MXl , MV , NAME , HAMEV , NFOUT) 
Ftinction: calculate frequency response from matrices 

General reference ?i2.3 

Response calculation: for last MX states only 


W 

A(MX2*MX2) 

ASQ(MX2*MX2) 

BO(MX*MV) 

MX 

MXl 

MV 

NAME(HX2) 

NAMBV(MV) 

NFOUT 


frequency (dimensionless) 
coefficient matrix A 

2 

coefficient matrix squared, A 
control matrix 

number of degrees of freedom 

number of first order degrees of freedom 

number of controls 

(maximum MX2= 2«MX - MX1= 60; maximiim MV = 6o) 
vector of variable names 
vector of control names 
unit number for printed output 
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STATIC 


Name i STATIC( A ,B0,MX,MX1,MV,NAME, NAMBV , NFOUT) 
Function t calculate static response from matrices 
General reference! section ?, 2,2 

Response calculation! for last MX states only 


A(MX2*MX2) 

bo(mx*mv) 

MX 

MXl 

MV 

NAME(MX2) 

NAMEV(MV) 

NFOUT 


coefficient matrix 
control matrix 

number of degrees of freedom 

number of first order degrees of freedom 

number of controls 

(maximum MX2= 2*MX - MXl = 60; maximum MV = 60) 
vector of vsuriable names 
vector of control names 
unit number for printed output 
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ZERO 


Name: ZER0(A ,B0,MX2,MX,MV ,NX,NV) 


Function! 

calculate zeros 


General reference: section 7.2.4 


A(MX2^^X2) 

coefficient matrix 


bo(mx*mv) 

control matrix 


MX2 

number of states f maximum 60 


MX 

number of degrees of freedom 


MV 

number of controls 


HX 

state number i for which zeros 

to be calculated 

NV 

control number j for which zeros to be calculated 

Output: 

LAMDAZ (MZ) 

zeros of x./v. 

^ j / iT 1 

[z-s) 

K1 

factor Kj^! 


MZ 

number of zeros 



available in the following common blocki 

GOMMON /EIGVZ/LAMDAZ(60) ,K1 ,MZ 
COMPLEX LAMDAZ 
REAL K1 
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ZETRAN 


Name: ZETRAN(Z,MZ) 

Function: transform matrix for zero calculation 

General reference: section 7.2.4 


Input: 

Z(MZ*MZ) matrix A* (A with column replaced by v. column of B) 

MZ number of states, MX2 


Output: 

z(mz*mz) 

MZ 


matrix A. (eigenvalues of which are the zeros); 
the facter is in Z(MZ*MZ+1) 

number of zeros 


GT 0 
EQ 0 
LT 0 


finite number of zeros exists 
no zeros, = Z(l) 
x^ not controllable by v. 
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BODE 


Name: B0DE(MX,MX1 ,MV,A2,A1,AO,BO,DOF1,NAMKX,NAMEV,NPLOT,NAMEXP,NAEMVP, 

NX , NV , NFO , NFl , ND.MSCALE , NFOUT) 

Function; calculate and printer-plot transfer function (Bode plot) 
General reference: section 7.2.3 


A2(MX*MX) 

Al(MX*Mxj 

A0(MX*MX) 

bo(mx*mv) 

MX 

MXl 

MV 

DOFl(MX) 

namex(mx) 

namev(mv) 

NPLOT 

NAMEXP(NX) 

NAMEVP(NV) 

NX 

NV 

NFO 

NFl 

ND 

MSCALE 

NFOUT 


coefficient matrices 


control matrix 

number of degrees of freedom 

number of first order degrees of freedom 

nvmiber of controls 

(maximum MX2= 2*MX - MXl = 60; maximum = 6o) 

integer vector designating first order degrees of 
freedom; DOFl(l) = 0 for x^^ first order 

vector of variable names 

vector of control names 

frequency response calculation method; if 1, from 
matrices; if 2, from poles and zeros 

vector of variable names to be plotted (inconsistent 
names ignored) 

vector of control names to be plotted (inconsistent 
names ignored) 

number of degrees of freedom to be plotted; meixlmum 30 

number of controls to be plotted; maximum 30 

exponent (base lO) of beginning frequency 

exponent (base lO) of end frequency 

frequency steps per decade 

(maximum NF= (NFl - NF0)*ND+ 1 = 151) 

magnitude plot scale; if 1, plot relative maximum value; 
if 2, plot relative 10**K; if 3 plot relative 10, 

unit number for printed output 
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BOISFP 


Name* BOIEPP(HM, HP, NFO.NFl.ND, OPTION, NFOUT) 

Functions printer-plot transfer function magnitude and phase 

HM(n) transfer function magnitude 

HP(n) transfer function phase (degrees, -180 to 180) 

(N = (NFl- NF0)*ND+1) 

NFO exponent (base lO) of beginning frequency 

NFl exponent (base lO) of end frequency 

ND frequency steps per decade 

OPTION magnitude plot scale: if 1, plot relative maximum 

value; if 2, plot relative 10**K; if 3» plot 
relative 10. 

NFOUT unit number for printed output 
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TRACKS 


Name i TRACKS( A 2 , A 1 , AO ,B0 ,MX , MXl ,HV , LOFl .CMEGA , NAMEX , NAMEV , NPLOT, 

PERICD, BELT, TMAX , NAMEXP.NAMEVP, NX , N\' , NFOUT) 

Fvinction: calculate and printer-plot time history of time- invariant 

system response 

General references section 7.2.5 

Calculate eigenvalue matrix and modal matrix: 

MRED = M without unused states (rows) 

f€B = without unused controls (columns) 

A2(MX*MX^ coefficient matrices 

Al(MX*MXj 

A0(MX*MX) 

B0(MV*MX) control matrix 

MX number of degrees of freedom 

MXl number of first order degrees of freedom 

MV number of controls 

(maximum MX2= 2*MX - MX1= 60; maximvun MV = 60) 

DOFI(HX) integer vector designating first order degrees of 

freedom; DCFl(l) = 0 for x^^ first order 

NAMEX(MX) vector of variable names 

NAfffiV(MV) vector of control names 

OMEGA frequency scale (rad/sec) 

NPLOT control input type 

1 step 

2 impulse 

3 cosine impulse 
k sine doublet 

5 square impulse 

6 square doublet 

PERIOD period T (sec) for impulse or doublet (NPLOT = 3 to 6) 

BELT time step (sec) 

TMAX maximum time (sec) 

(maximum NX*NV*TMAX/BELT = ?200) 
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TRACKS 


namexp(nx) 

NAMEVP(NV) 

NX 

NV 

NFCUT 


vector of variable names to be plotted (inconsistent 
names ignored) 

vector of control naunes to be plotted (inconsistent 
names ignored) 

number of degrees of freedom to be plotted; maximum 30 
number of controls to be plotted; maximum 30 
unit niimber for printed output 
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TRCKPP 


Name ; TRC3CPP( TRA CE , NX , NV , MT , CELT , NAMEXP , NAMEVP , NFOUT ) 
Functio'': printer- plot time history 


TRACE(NX,NV,MT) 

NX 

NV 

MT 

BELT 

NAMEXP(NX) 

NAMEVP(NV) 

NFOUT 


array of time history traces to he plotted 
number of degrees of freedom to he plotted 
number of controls to he plotted 
(maximum NX*NV = 26) 
number of time steps to be plotted 
time step (sec) 
vector of variable names 
vector of control names 
unit number for printed output 


i 
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GUSTS 


Name: GUSTS(A2,A1,AO,BO,MX,MX1,MV,MG,DOF1,NAHEX,RADIUS,OMEGA,GRAV, 

EULER , VEL , LGUST , MGUST , HA ME® , HAMEXL , ML , NAMEX A , MA CC , 
FREQA ,RACC , HEM , ZETA , NAMEXB ,NFCUT) 

Functionj calculate and print rms gust response 

General reference: section 7.2.6 


A2(MX*MX^ 

A1(MX*MX) 

AO(MX*MX) 

BO(MX*M\0 

MX 

MXl 

MV 

MG 


K)Fl(MX) 

NAMEX(MX) 

RADIUS 

OMEGA 

(®AV 

EULER(2) 

VEL(3) 

lgust(mg) 


mgust(mg) 

NAMEXR(3) 


coefficient matrices 


control matrix (gust in last MG columns) 

number of degrees of freedom 

number of first order degrees of freedom 

number of controls and gusts 

number of gust components 

(maximum MX2 = 2*MX - MXl + MAGC + MG = 60) 

( maxim lom MG = 3) 

integer vector designating first order degrees 
of freedom; DOFl(l) = 0 for first order 

vector of variable names 

length scale R (ft or m) 

frequency scale iX (rad/sec) 

acceleration due to gravity (ft/sec or m/sec ) 

trim Euler angles and (rad); required 

for body axis acceleration only 

velocity components in body axis frame (divided by 
IXR); only magnitude required (for *®q) unless body 
axis acceleration calculated 


real vector of gust correlation lengths: if GT 0, 

dimensional correlation length L ( 'Cq = L/2V); if 
EQ 0, L = 400. used; if LT 0, magnitude is correlation 
time (dimensionless), so break frequency is 

real vector of gust component relative magnitudes 


Ic 


in state vector (NAMEX); 


names of , 

analysis assumes that t*1s’ ^Is’ '^Is 

immediately (inconsistent names ignored) 


follow 
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GUSTS 


namexl(ml) 

ML 

nambxa(macc) 

freqa(macc) 

MACG 

RACC(3) 

ZETA(3,NEM) 

NEM 

NAMEXB(6f+NEM) 

NFOUT 


names of linear degrees of freedom in state 
vector (NAMEX) for dimensional output (ft or m, 
obtained from R)i degrees of freedom not 
identified are angular (degrees) (inconsistent 
names ignored) 

number of linear degrees of freedom 

names of degrees of freedom ( NAMEX ) for which 
acceleration calculated; last three names must 
eaual ACCB to calculate body axis acceleration 
(all three or none) (inconsistent names ignoredj 

accelerometer break frequency (Hz), in same order 
as NAMEXA; 2/rev used if FREQA ^ 0. 

number of accelerometers; none if MACC $ 0 

X, y, z location of point at which bodv axis 
acceleration calculated (dimensionless) 

airframe elast mode shapes, k - 1 to NEM; required 
for body axis acceleration only 

number of airframe elastic modes; none if NEM ^ 0; 
maximum 10 

names of S^p. *p, yp, Zp, qp^ ••• ‘IFnem 

in state vector (NAMEX); assumes all elastic airframe 
states are consequtlve; required for body axis 
acceleration only (inconsistent names ignored; 

unit number for printed output 
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PSYSAN 


Namei PSYSAN( MX, MX1,A2,A1, AO, PHI, DT, NT, MT, PERIOD, DOFl, HINT, NFOUT) 

Function: analyze system of periodic coefficient linear differential 

equations 

General reference: section 7.3 


A2(MX*MX) 

ai(mx*mx) 

A0(MX*MX) 

MX 

MXl 


DOFl(MX) 


DT 

NT 

MT 

PEKIOD 

PHI 


NINT 


NFOUT 


coefficient matrices 


number of degrees of freedom 

munber of first order degrees of freedom 

(maximum MX2 = 2*MX - MXl = 60) 

integer vector designating first order degrees 
of freedom (zero colximns in AO); DOFl(l) = 0 
for Xj^ first order 

time increment; may vary with NT, but for Runge-Kutta 
integration successive pairs must be eq\al 

time step counter (NT =0, 1, 2, ... MT) 

total nianber of time steps in numerical integration; 
for Runge-Kutta integeration, must be even 

period T of the system 

temi»rary storage of state transition matrix ^ and 
last A; dimension 2*MX2*MX2 for modified trapezoidal 
integration; dimension 3*MX2*MX2 for Runge-Kutta 
integration (MX2» 2*MX - MXl) 

numerical integration method: if 1, modified 

trapezoidal method, error order DT**3» if 2, 
Runge-Kutta method, error order (2*LT)**5 

unit number for printed output 


Output: 

LAMDA(MX2) 

LAMDAC(MX2) 

MX2 


roots V (principal value) 
eigenvalues of i(T) 
number of poles 

available in the following common block: 

COMMON /EIGVP/UMDA(60) ,LAMDAC(60) ,MX2 
COMPLEX LAMDA,LAMDAC 
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PSYSAN 


Typical usage: 

DT =- PERIOD/MT 
DO 1 NT “ O.MT 
T - DT * NT 

calculate coefficient matrices at time T 
1 GALL PSYSAN 
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EEIPRAN 


Nanet EEPRAN(A,MX,MX1,A2,A1,A0,D0F1,NF0UT) 

Function: transform equations to state variable form 

General reference: section ?.l 

Input: 

A2(MX*MX 
A1(MX*MX 
AO(MX^TOC 

MX 

MXl 

DOFl(MX) 

NFOUT 

Output: 

A(MX2*MX2) coefficient matrix (MX2 = 2*MX- MXl) 


coefficient matrices 

number of degrees of freedom; maximum 60 

number of first order degrees of freedom 

Integer designator of first order degrees 
of freedom; IDFl(l) “ 0 for first order 

unit number for printed output 
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MAINTB 


Name I MAINTB 

Function* airfoil table preparation 
General reference* section 2«4.4 

Subprograms required* AEROTi AEROPPi C81INT« C8lRD» REDCL, TABFIX 
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AEROT 


Namei AER0T(ALPHA, MACH, RADIAL, OPTION, a, CD, CM) 
Function: interpolate airfoil tables 

General reference: section 2,4,4 


ALPHA 

MACH 

RADIAL 

OPTION 

CL 

CD 

CM 


angle of attack (deg) 

Mach number M 
radial station r/R 

integer parauneter: if 1 calculate Cj| j if 2 

calculate c , , if 3 calculate c , if 4 calculate 
all three coefficients 

°d2D 

°®2D 
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AEROPP 


Name: AEROPP(CL,GD,GM,MA,AMAX) 

Function: printer-plot airfoil aerodynamic cheuracteristics 

Galculate ordinate limits: 



a) 

c = maximum value of magnitude 


b) 

N= [log cl (N = N- 1 if c-< 1.) 


c) 

K = [,c/10**n1 + 1 


ci) 

use for scale X = K * 10**N 

cl(ma) 


array of to be plotted 

gd(ma) 


array of to be plotted 

gm(ma) 


array of c to be plotted 
m 

MA 


number of angle of attadc values; odd number 

AMAX 


maximum angle of attack; data in arrays for 

eu = — to , in MA steps 

max max’ 
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3. CX)MPinER SYSTEM SUBPROGRAMS 


The following computer system subprograms (or the equivalent) 
are required to determine the calendar date and time of day, which 
form the identification for jobs and files. 

a) CALL TIME(ITIME) 

Function; returns time of day (8 alphanumeric 
characters) in array ITIME(2) 

b) CALL DATECIDAIE) 

Function: returns calendar date (8 alphanumeric 

characters) in array IDA'IE(2) 

The following computer system subixrograms (or the equivalent) axe 
required in the timing subprogram. 

a) CALL SETTIM(O.O) 

Function: initializes timer 

b) ITIME = INTVAL(O.O) 

Function: returns CPU time, in milliseconds since 

initialization 



4 CORE REQUIREMENTS 


The program requires 4.04 megabytes of core storage. Of this 
total, 1.84 megabytes is for the subprograms and 2.20 megabytes is 
for the common blocks. The common blocks for the nonuniform inflow 
influence coefficients (both rotors) require 0.96 megabytes. 
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